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Bi-spectral  satellite  thresholds  for  precipitation  speci¬ 
fication  are  explored  with  visual  and  infrared  satellite  data 
collocated  with  Service-A  hourly  observations  for  137  surface 
stations  in  the  southeastern  United  States.  The  data  span  the 
month  of  August  1979  and  total  70,623  observations,  including 
538  daylight  precipitation  observations. 
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explored  and  determined  to  be  significant  in  the  representation 
of  weather  conditions.  Precipitation  and  no-precipitation  data 
can  be  statistically  differentiated  with  the  visual  and  infrared 
mean  and  standard  deviation  values. 

For  overcast  ceiling  reports,  a  simple  linear  bi-spectral 
threshold  based  on  a  50%  probability  of  precipitation  is  defined 
as  extending  from  albedo  1.00  to  0.60  with  associated  cloud  top 
temperatures  290K  and  210K,  respectively.  For  overcast  and  bro¬ 
ken  ceiling  reports,  an  albedo  greater  than  0.80  specifies  a  50% 
probability  of  precipitation. 
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ABSTRACT 

-spectral  satellite  thresholds  for  precipitation  spec¬ 
ification  are  explored  with  visual  and  infrared  satellite 
data  collocated  with  Service-A  hourly  observations  for  137 
surface  stations  in  the  southeastern  United  States.  The 
data  span  the  sonth  of  August  19  79  and  total  70,623  observa¬ 
tions,  including  538  daylight  precipitation  observations. 

The  distributional  and  statistical  differences  of  four 
satellite  resolution  sizes  ranging  fros  484  to  2025  nai*  are 
explored  and  deterained  to  be  significant  in  the  representa¬ 
tion  of  weather  conditions.  Precipitation  and  no-precipita- 
tion  data  can  be  statistically  differentiated  with  the 
vi-sual  and  infrared  aean  and  standard  deviation  values. 

For  overcast  ceiling  reports,  a  siaple  linear  bi-spec- 
tral  threshold  based  on  a  50X  probability  of  precipitation 
is  defined  as  extending  froa  albedo  1.00  to  0.60  with  asso¬ 
ciated  cloud  tcp  teaperatures  2  90K  and  21  OK,  respectively. 
For  overcast  and  broken  ceiling  reports,  an  albedo  greater 
than  0.80  specifies  a  50%  probability  of  precipitation. 
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I.  INT  ROD  OCT  ION 


The  determination  of  precipitation  occurrence  and 
amounts  is  an  important  factor  in  scientific,  commercial, 
and  operational  endeavors.  Scientific  uses  are  concentrated 
in  the  fields  cf  meteorology,  hydrology,  and  oceanography, 
where  precipitation  is  essential  in  analysis,  diagnosis, 
prediction,  and  verification.  Hithin  meteorology,  precipi¬ 
tation  serves  as  both  a  forcing  and  response  element  in  the 
study  of  daily  weather  and  climatology.  Indeed,  precipita¬ 
tion  is  a  critical  input  for  climate  research  and  into  gen¬ 
eral  circulation  models  which  premise  to  extend  the  time 
frame  of  skillful  weather  forecasts.  Commercial  uses  encom¬ 
pass  agriculture,  forestry,  transportation,  communications, 
water  resource  management,  and  many  others. 

Despite  the  importance  of  precipitation  data  to  a  vari¬ 
ety  of  fields,  there  are  serious  shortcomings  in  current 
precipitation  determination.  These  shortcomings  are  due  to 
areal  and  economic  limitations  imposed  upon  the  land-based 
rainfall  monitoring  systems.  h  possible  solution  is  embod¬ 
ied  in  precipitation  information  extracted  from  satellite 
data.  With  the  advent  cf  high  resolution,  multi-spectral 
channel  satellites  in  the  late  1970 's,  satellite  derived 
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precipitation  data  are  being  studied  as  a  viable  method  to 
complement  and  supplement  conventional  rainfall  data. 

The  satellite  image  interpreter  does  a  subjective  analy¬ 
sis  based  on  the  gray  shade  variations,  representing  a  range 
of  digital  counts,  that  appear  in  the  satellite  image.  How¬ 
ever,  satellite  data  contain  more  information  within  the 
digital  values  than  can  be  resolved  by  the  human  eye  in  pho¬ 
tographic  images.  The  satellite  digital  counts  input  into  a 
computer  allow  use  of  the  full  range  of  the  digital  values. 

Until  recently,  computer  processing  of  satellite  data 
has  been  confined  largely  to  research  uses.  Acquisition, 
storage,  and  processing  of  the  huge  volumes  of  digital  sat¬ 
ellite  data  could  not  be  handled  operationally  in  real  time. 
However,  with  the  recent  advent  of  more  capable  mini-com¬ 
puter  systems,  such  as  the  United  States  Navy's  Satellite 
Data  Processing  And  Display  System  (SPADS)  developed  by  the 
Naval  Environmental  Prediction  Research  Facility  (NEPRF) 
Honterey,  California,  real  time  quantitative  use  of  digital 
satellite  data  has  become  a  reality.  Bith  the  operational 
availability  of  such  systems  as  the  SPADS  unit,  there  is  a 
need  for  numerical  schemes  to  aid  in  the  objective 
specification  of  current  weather  conditions. 
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This  thesis  concentrates  on  the  specification  of  visual 
and  infrared  satellite  data  thresholds  in  etereiniag  pre¬ 
cipitation  occurrence  and  qualitative  precipitation  rates  in 
a  aid-latitude  coastal  environaent.  The  data  set  used  con¬ 
sists  of  collocated  Geostationary  Operational  Environaental 
Satellite--East  (GOES-E)  satellite  data  and  hourly  surface 
observations  at  East  coast  and  Gulf  Coast  United  states  sta- 
tions,  south  of  40°l,  for  the  aonth  of  august  1979. 

The  use  of  satellite  data  for  precipitation  specifica¬ 
tion  is  not  new.  There  is  the  recognized  linitation  that 
infrared  and  visual  satellite  sensors  are  neasuring  proper¬ 
ties  associated  with  snail  cloud  particles  and  not  precipi¬ 
tation  sized  particles.  fonetheless*  Buench  and  Keegan 
(1979)  specified  quantitative  precipitation  rates*  Liljas 
(1981a*  1981b)  specified  qualitative  precipitation  rates* 

and  love  joy  and  Austin  (1979)  delineated  rain  versus  no-rain 
cases  using  visual  and  infrared  satellite  data.  Del  Beato 
(1981)  used  cloud  top  tesperatur es  within  a  restricted  cloud 
case  classification  '  derive  qualitative  precipitation 
rates. 

The  currently  avai.1^  <=  precipitation  study  results  are 
based  on  data  sets  vith  region*  season*  and  size 


limitations.  This  research  effort  will  use  data  from 
stations  covering  aore  than  42  0,000  square  nautical  ailes 
(nmi*)  in  the  eastern  and  central  United  States  with  a  total 
of  70,623  observations  (538  precipitation  observations)  .  In 
coaparison,  the  relatively  comprehensive  precipitation  study 
of  Nuench  and  Keegan  (1979)  was  based  on  552  cases  (300 
rainfall  cases)  froa  five  stations  in  the  northeastern 
United  States  for  April  through  November  1977.  The  signifi¬ 
cantly  larger  size  of  the  present  sample  will  allow  better 
statistical  determination  of  appropriate  distributions  and 
threshold  value  significance. 

The  primary  objective  of  this  thesis  is  to  investigate 
specification  of  precipitation  versus  no- precipitation  froa 
satellite  visible  and  infrared  digital  counts,  Addition¬ 
ally,  in  precipitation  cases,  the  feasibility  of  qualitative 
specification  of  light  versus  aoderate/heavy  precipitation 
and  quantiative  specification  of  convective  versus  continu¬ 
ous  precipitation  are  investigated. 

Th€  thesis  is  organized  into  five  chapters.  Chapter  II 
reviews  the  satellite  data  based  precipitation  studies. 
Chapter  III  describes  the  data  set,  the  data  processing  and 
the  testing  program.  Chapter  IV  describes  the  results. 
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Chapter  7  states  the  conclusions  and  suggests  further 
research. 


II.  PRECIPITATION  SPECIFICATION 


A.  INTRODUCTION 

The  specification  of  precipitation  and  the  estination  of 
rainfall  rates  using  satellite  isagery  have  bean  studied 
using  a  variety  of  aethods  over  a  vide  spectrua  of  tine 
scales.  This  review  vill  concentrate  on  those  methods 
developed  for  synoptic  scale  and  nesoscale  analysis  of  pre¬ 
cipitation  on  a  diurnal  or  shorter  tine  scale.  The  aethods 
reviewed  include  bi-spectral  and  infrarei  threshold  (Huench 
and  Keegan,  1979;  Liljas,  1981a,  1981b;  Lovejoy  and  Aus¬ 
tin,  1979;  Del  Beato,  1981;  Rylie,  1982)  and  life  history 
(Scofield,  1981;  Griffith  et  al. ,  1978;  Stout  jst  al. , 

1979;  Uyli9,  1979;  Negri  and  Adler,  1981). 

B.  BI-SPECTRAL  AND  INFRARED  THRESHOLD 

The  bi-spectral  threshold  method,  in  which  infrared  and 
visual  satellite  data  are  used,  involves  napping  the  extent 
and  distribution  of  precipitation.  Coabining  the  visual  and 
infrared  data  provides  information  on  the  cloud  temperatures 
(infrared  data)  and  on  the  cloud  thickness  (visual  data). 
Thus,  while  use  of  the  visual  or  infrared  data  alone  may 


have  limitations  in  specifying  precipitation,  the  coabina- 
tion  of  both  aay  succeed  at  specifying  precipitation.  The 
aulti-spectral  satellite  channels  introduced  on  satellites 
in  the  late  1970*s  yielded  the  possibility  of  bi-spectral 


thresholds . 

Threshold 

values  and  study 

condition 

paraaeters 

of  selected 

bi-spectral 

studies 

are  suaaarized  in 

Table  I. 

TABLE 

I 

Suaaary 

of  Bi-spectral  and  Infrared 

Threshold 

Values 

STOW 

LOCATION 

CASES 

TIME  OF 
TEAK 

THRESHOLD 

INFRARED 

VALUE 

VISUAL 

(19TO) 

1  Northeastern 
United  Ststns 

532 

obser¬ 

vations 

April- 

lonabir 

1977 

-12°C 

0.60 

Llljaa 

(1981) 

Scsndlasrls 

- 

May  1979 
August 
1979 

-12°C  to 
-15°C 

- 

Lora  Joy  and  Hoacraal 

Austin 
(1979) 

17  days 

Jaaa  1977 

-21°C. 

-26°C. 

-41°C 

.80.  .88* 

*  Visual  threshold  baaod  on  norasllxsd  acala  from  0-1 

Huench  and  Keegan  (1979)  studied  precipitation  specifi¬ 
cation  using  GOES  visual  and  infrared  satellite  data  and 
hourly  rainfall  cliaatological  data  for  five  stations  in  the 
northeastern  Onited  States  for  the  period  April  through 
Hoveaber  1977.  Their  data  set  consisted  of  552 
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observations,  coaprised  of  300  rainfall  observations  and  252 
cases  of  either  aonprecipitating  cloudy  or  fair  weather 
observations.  The  visual  (1  ka  resolution)  and  infrared  (7 
ka  resolution)  satellite  data  were  area  averaged  over  7x7 
square  kiloaeters  (ka2)  and  14  x  14  ka2,  respectively.  A  65 
point  visual  data  array  (8x8  plus  the  center  point)  and  a 
17  point  infrared  data  array  (4x4  plus  the  center  point) 
were  centered  ever  each  station.  The  GOES  visual  data  were 
noraalized  using  reflection  values  froa  Liou  (1976)  with  the 
aodification  of  lower  absorption  and  higher  transaission  to 
coapensate  for  Liou's  treataent  of  the  couplets  solar  spec** 
trua.  The  anisotropic  radiation  of  clouds  was  corrected 
with  functions  calculated  by  fiuench  and  Keegan  (1979)  from 
ground-based  radiometers  and  satellite  measurements.  Froa 
these  data,  they  determined  probabilities  for  precipitation 
greater  than  .01  and  .10  inches  in  one  hour  and  the  amount 
of  precipitation  for  the  hour  following  the  satellite 
observation  (see  Figs.  1,  2,  and  3). 

Huench  and  Keegan  (1979)  did  not  provide  the  standard 
deviations  for  the  data  in  these  figures.  However,  they 
stated  there  was  "considerable  uncertainty  in  the  specifica¬ 
tion  of  rainfall  amount."  As  an  example,  they  stated  that 
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Figure  3.  Probability  of  One-hour  Rainfall  Greater  than 

0. 10  Inches  as  a  Function  of  Cloud  Reflectivity 
and  Infrared  Cloud  Temperatura  (from  Huench  ana 
Keegan*  1979) 

using  Fig.  3  "for  a  one-hcur  rainfall  specification  of  0.  10* 
two- thirds  of  the  values  would  fall  between  0.25  and  0.0  4." 
Huench  and  Keegan  stated  that  their  figures  emphasize  the 
reguireaent  fcr  both  visual  and  infrared  data  to  specify 
precipitation  asounts. 

Liljas  (1981a*  1981b)  developed  a  bi-spectral  cloud 

classification  based  on  visual  and  infrared  data  from  the 
polar  orbiting  TIB0S--6  satellite  (see  Fig.  4  and  Table  II)  . 
The  data  set  consisted  of  a  limited  number  of  daily  observa¬ 
tions*  chosen  for  their  synoptic  characteristics*  in  Hay  and 
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Figure  4.  Two  Dimension  Decision  Space  for  Typing  Clouds 

free  Visual  and  Infrared  Digital  Counts  (Table  II 
defines  the  symbols  for  the  clouds)  (from  Liljas, 
1981a) 


August  1979  over  a  region  encompassing  Norway,  Sweden,  Fin¬ 
land,  and  the  Baltic  Sea  with  weather  charts  providing  the 
ground  truth.  Based  upon  the  precipitation  threshold 
results  of  Buench  and  Keegan  (1979),  Liljas  chose  a  cloud 
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TABLE  II 

Cloud  Classification  to  be  used  with  Fig.  4 
(froa  Liljas,  1981a) 


Mala  and  Cloud  Typeai 

1.  Cumulonimbus 

2.  liaboatracua 

3.  Clrroacraeua 


4.  Cuaulua  eongaaeaa 


3.  Stracociwlu* 


4.  Baze/Scr atua 


7. 


I.  Bacas 


b 

* 

e 

t 

t 

1 

h 

I 

k 

l 

Z 

1 

1 


Scon  cloud  with  high  top 

Squall  cloud  vlch  scattered  ahowara 

Largo  ware leal  chlckaaaa 
lather  low  eopslda 

Daaaa  clrroacraeua 
Clxrua 

Thin  elrrua  over  wacar 
Daaaa  alcoacracua 
Thin  alcoacracua  ovor  wacar 
Thin  alcoacracua 

Daaaa  alcoeuaulua 
Large  piled  up  cuaulua 
lacher  nail  and  £lac  cauloa 


Daaaa  scratocuaulua 
Ordinary 

Slightly  piled  up  cuaulua 
with  dear  araaa  in  between 
▼ary  dense  haze/scracua 
Daaaa  haze/scracua 

Ordinary  haze/scracua 
Cinulis  humilis 


Baza  over  water 


1.  Planting  season  spring 
1^  Warn  green  season 


or  autumn 


k 

k.  Cold 
kj  Vara 
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top  temperature  threshold  of  -,12°C  to  -15°c  to  classify 
cumulonimbus  and  niabostratus  clouds.  starting  with  this 
cloud  classification  and  the  assuaption  that  the  highest  and 
densest  clouds  produce  the  maximum  precipitation  amount, 
Liljas  suggested  a  qualitative  precipitation  intensity  scale 
based  on  the  sub  of  the  visual  and  infrared  satellite  digi¬ 
tal  counts  (see  Table  III)  .  These  sums  represent  the  areas 
of  the  Liljas  niabostratus  and  cumulonimbus  cloud  types  in 
his  bi-spectral  cloud  classification  (see  Pig.  5). 

TABLE  III 

Threshold  Values  Describing  Precipitation  Intensity  Levels 
as  Applied  in  Pig.  5  (from  Liljas,  1981a) 

The  Sum  al  Digital  Levels 

Ch  1  ♦  Ch  4  391-310  light  rain 

31 1-330 
331-350 
35  1-370 
371-390 

390  very  strong  rain 

Lcvejoy  and  Austin  (1979)  studied  rain  mapping  of  cloud 
areas  based  on  GOES  visual  and  infrared  satellite  data  over 
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Figure  5.  Precipitation  Intensi  ty ,  Classification  from 
Visual  and  Infrared  Digital  Counts.  The 


mathematical  descrip 
(from  Liljas,  1981af 


Montreal ,  Canada,  and  the  tropical  Atlantic  (Global  Atmos¬ 
pheric  Research  Program  Atlantic  Tropical  Experiment,  GATE, 
data)  with  radar  data  providing  the  ground  truth.  The  Mont¬ 
real  data  set  consisted  of  17  observations  over  three  days 


daring  June  1977.  Working  with  4  x  4  km  resolution  satel¬ 
lite  iaage,  lovejoy  and  Austin  plotted  two  dimensional  fre¬ 
quency  grids  for  the  radar-de termined  rain  and  no-rain 
points  on  a  25  x  25  array  (see  Figs.  6  and  7).  The  visual 
data  were  normalized  by  selecting  the  ••brightest"  and  "die- 
nest"  values  in  each  iaage  and  linearly  interpolating  the 
radiances  between  0  and  1. 

Lovejoy  and  Austin  (1979)  state,  with  reference  to  the 
cunulus  rain  data  distribution  of  Fig.  6  that,  "The  distri¬ 
bution  was  to  a  good  approximation  a  two-diaensional  Gaus¬ 
sian."  They  do  not  describe  or  provide  the  statistics  to 
support  this  assertion.  The  no-rain  cumulus  cases  (Fig.  7) 
were  described  as  a  biaodal  distribution  with  one  peak  near 
the  low  visual  and  low  infrared  values  and  the  other  peak 
near  the  rain  peak  but  shifted  slightly  toward  lower  values. 
In  most  cases,  the  separation  of  the  oumulus  rain  and  no¬ 
rain  cases  was  statistically  significant  with  the  probabil¬ 
ity  ranging  from  10 *  to  50%  that  the  rain  and  no-rain 
samples  came  from  the  same  population. 

The  Lovejoy  and  Austin  (1979)  two  dimensional  frequency 
plots  for  non-cuaulus  storms  were  limited  to  one  case.  The 
significant  differences  between  the  cumulus  and  non-cumulus 
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Figure  6. 


Frequency  Plot  of  Bain  Distribution  for  GATE  day 
2  4  8,  1300  GST  (from  Lovejoy  and  Austin,  1979) 
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Figure  7. 


Frequency  Plot  of  No-rain  Distribution  for  GATE 
day  248,  1  300  GST  (from  Lovejoy  and  Austin,  1979) 
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data  sets  war*  that  the  non-cumulus  no-rain  plot  lost  its 
bieodal  character,  relative  to  the  cumulus  no-rain  plot,  and 
appeared  as  a  broad  two  diaensional  Gaussian  distribution. 
The  non-cuaulus  rain  plot  points  fell  within  the  no-rain 
distribution,  but  were  shifted  slightly  higher  in  the  vis¬ 
ual.  The  separation  of  the  non-cuaulus  rain  and  no-rain 
cases  was  not  statistically  significant,  with  greater  than  a 
SOX  probability  of  the  rain  and  no-rain  saaples  coming  from 
the  same  population. 

Love  joy  and  Austin  (1S79)  attempted  to  further  classify 
the  cumulus  rain  and  no-rain  cases  into  no-rain,  light  rain, 
and  heavy  rain.  Rainfall  rates  greatar  than  2  mm-h-*1,  as 
determined  by  radar,  were  defined  as  heavy  rain.  As 
expected,  the  mean  of  the  heavy  rain  cases  was  shifted 
slightly  towards  higher  visual  and  infrared  values  than  the 
mean  of  the  light  rain  cases.  However,  the  shift  was  so 
small  that  there  was  at  least  an  80%  probability  of  the 
light  rain  and  heavy  rain  cases  coming  from  the  same  popula¬ 
tion.  Lovejoy  and  Austin  (1979)  concluded  that  "little  if 
any  rainfall-rate  information  is  contained  in  a  single  (vis¬ 
ual  and  infrared)  satellite  image." 


Love joy  and  Austin  (1979)  tested  a  spectral  threshold 
technique  for  rain  area  napping.  Each  satellite  image  of 
400  z  400  ka  was  divided  into  one  hundred  40  z  40  km  bozes. 
The  100  subareas  were  each  checked  with  radar  to  determine 
the  total  number  of  rain  areas.  An  equal  total  number  of 
satellite  subareas  were  classified  as  rain  areas.  The  sat¬ 
ellite  subareas  with  the  highest  visual  and  highest  (cold) 
infrared  values  ware  classified  as  rain  areas,  until  the 
total  number  cf  satellite  rain  areas  equaled  the  total  num¬ 
ber  of  radar  determined  rain  areas.  This  spectral  threshold 
technique  was  applied  to  three  days  accumulation  of  data  and 
is  shown  in  Tables  IT  and  V.  When  compared  with  the  success 
of  the  two  dimensional  frequency  plot  method,  the  visible 
and  infrared  thresholds  averaged  45%  and  58%  worse,  respec¬ 
tively.  The  accuracy  of  the  visual  threshold  is  limited  by 
the  eztent  of  low,  thick  clouds  and  the  infrared  threshold 
is  limited  by  the  eztent  of  the  cirrus  clouds  in  the  satel¬ 
lite  image.  Lovejoy  and  Austin  (1979)  concluded  that  "the 
errors  involved  in  using  a  'best  threshold'  are  very  large 
indeed. " 

Del  Beato  (1981)  studied  correlations  between  cloud  top 
temperatures  (based  on  HOAA-5  satellite  data)  and  rainfall 


TABLE  IV 

Statistical  Comparison  of  Rain  Area  Sapping  Techniques 


(R  /R  X  100  indicates  "percentage  of  correct  satellite  rain") 
(froe  Lovejoy  and  Austin,  1979) 
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TABLE  V 

Statistical  Comparison  of  the  Accuracy  of  Rain  Areas 
(from  Lcvejcy  and  Austin,  19  79) 


Technique 

Region 

Number  of  Images 
or  Sequences 

Bias 

Error 

Factor 

E.M. 

2-D  Pattern 
Matching 

Montreal 

17 

1.13 

1.26 

0.22 

2-D  Pattern 
Matching 

Montreal 

3 

1.08 

1.19 

0.18 

Optimum  IR 
Threshold 

Montreal 

3 

1.38 

1.74 

0.71 

Optimum  Visible 
Threshold 

Montreal 

3 

1.54 

1.59 

0.58 

2-D  Pattern 
Matching 

CATE 

S 

1.21 

1.41 

0.25 

totals  for  30-  and  60-oin  intervals  over  eastern  Australia. 
The  satellite  data  had  a  60  km2  maximum  resolution  at  subsa¬ 
tellite  point  and  cloud  top  temperatures  were  area  averaged 
for  a  resolution  of  200  km2.  The  21  data  sets  were  first 
classified  according  to  synoptic  situation  in  a  rough 
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attempt  tc  group  the  data  by  cl  cud  type,  droplet  spectra, 
and  air  ease  trajectory.  The  initial  results  suggested  that 
the  cloud  top  teaperature  determined  an  upper  limit  on  rain* 
fall  amount,  with  the  maximum  increasing  as  the  cloud  top 
teaperature  decreases.  k  linear  correlation  analysis  to 
determine  a  quantitative  relationship  between  rainfall 
amount  and  cloud  top  teaperature  gave  indefinite  results. 

Further  study  of  surface  and  radiosonde  observations 
indicated  that  classif ication  by  proportion  of  cumuliform 
cloud  reports  to  all  cloud  reports  and  subcloud  layer  humid¬ 
ity  might  be  acre  appropriate  (Del  Beato,  1981).  This  clas¬ 
sification  resulted  'in  a  correlation  coefficient  of  0.90, 
excluding  cases  with  cumuliform  portions  less  than  50%  and 
dew-point  depressions  of  greater  than  6°C.  Finally,  a  com¬ 
posite  frequency  distribution  was  calculated  based  on  three 
cases,  all  southwesterly  stream  situations  described  as 
"post-frontal  cellular  convection  cases  in  cyclonically 
curved  flow."  The  fitted  equation  was: 

f  *  0.C57  -  0.004CTT  -  0.054H  (1) 

where  f  is  the  rainfall  frequency,  B  is  the  30-min  rain 
total  (mm)  ,  and  CTT  is  the  cloud  top  temperature  (°C)  .  The 
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agnation  vas  fitted  to  41  independent  f  values.  This  equa¬ 
tion  is  associated  with  a  correlation  coefficient  of  0.79  at 
the  99%  confidence  level.  Equation  (1)  indicates  no  rain 
fros  clouds  varner  than  +13°C  and  a  maximum  30-ein  rainfall 
of  2.5  an  for  a  cloud  top  teaperature  of  -20°c. 

In  suaaary*  Del  Beato  (1981)  found  that  cloud  top  temp- 
eratures  and  30-  and  60-ain  rainfall  totals  indicated  sta¬ 
tistically  significant  relationships  for  cloud  systems  vith 
a  high  proportion  of  cuaulus  clouds  and  high  subcloud  humid¬ 
ity.  Additionally*  as  cloud  top  temperatures  decrease  to  at 
least  -35°C*  rainfall  totals  increase. 

Nylie  (1982)  atteapted  to  correlate  rainfall  occurrence 
vith  radiosonde  soundings*  hourly  Service-A  observations* 
and  visual  and  infrared  satellite  data.  His  data  sample  was 
restricted  to  "large-scale  cloud  cover"  areas  vith  wide¬ 
spread  precipitation  (rain  gauge  reports  varied  less  than 
20%)  for  the  Great  Plains  States  region  for  the  period  27 
February  1981  through  4  January  1982.  Froa  thirteen  parame¬ 
ters  derived  from  the  three  data  sources  (see  Table  VI) *  the 
best  linear  regression  equation  for  estimating  rainfall 
rates  vas: 

6-hour  rain  (in)  *  1.0242  «■  0.380PW  -  0.Q304QC 

-  0.00 47Ct  (2) 
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where  Pw  is  the  vertically  integrated  precipitable  water 
vapor  (in)  ,  Qc  is  the  moisture  convergence  (g/kg/day)  ,  and 
Ct  is  the  cloud  top  temperature  (Kelvins)  .  Equation  (2)  has 
a  linear  correlation  coefficient  of  0.60.  Linear  regression 
equations  were  also  deterained  for  the  three  parameters 
alone  and  for  a  combination  of  Pw  and  Qc  to  be  used  when  not 
all  three  data  types  were  available.  The  cloud  temperature 
regression  equation  was: 

6  hour  rain  (in)  *  2.10  -  0  .008Ct  (3) 

The  correlation  coefficient  was  -0.35.  Hylie  (1982)  stated 
that  the  synoptic  scale  data  base  measurements  were  best 
suited  for  estimating  broad  changes  in  rainfall  rates  asso¬ 
ciated  with  changes  in  air  masses  and  not  suited  for  esti¬ 
mating  rainfall  rates  associated  with  small  scale  dynamic 
processes. 

C.  LIPE  HISTORY 

The  life  history  methods  are  empirically  derived  precip¬ 
itation  estimation  schemes  based  upon  two  assumptions, 
first,  that  significant  rainfall  comes  from  convective 
clouds,  and  second,  that  convective  clouds  can  be  identified 
and  measured  in  satellite  images.  These  methods  involve 
manual  analyses  of  convective  cloud  areas  in  a  sequence  of 
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TABLE  71 

Correlation  Coefficients  for  Determination  of  Precipitation 


Based  on  the  Three  Data  Types  (from  Wylie,  1982) 


MEASURED 

E4MSJ2SH 

CORRELATION 

WITH  6  H00R 

MKI£.  beport 

NUMBER 

OF  CASES 

F 

Vertically  integrated 
precipitable  water  vapor 

0.48 

196 

5  8* 

Cloud  top  brightness 

-0.44 

184 

44* 

Cloud  top  height 

-0.40 

190 

36* 

Hoisture  convergence 

0.38 

184 

31* 

Cloud  top  temperature 

-0.35 

199 

27* 

Bubble  model  predicted  cond. 

0.27 

115 

9 

500  mb  vorticity  advection 

-0.21 

173 

8* 

Parcel  lifted  index 

-0.20 

200 

8* 

700  mb  temperature  advection 

0.20 

173 

7* 

Sfc  temperature  advection 

0.19 

156 

6 

850  mb  temperature  advection 

0.17 

189 

6 

Wind  convergence  (sfc) 

0.09 

167 

1 

Vertical  wind  shear 

0.03 

156 

0 

*  Significant  correlation  at  the  99X  level. 

visual,  infrared,  or  both  visual  and  infrared  satellite 
iaages.  Threshold  values  and  study  condition  parameters 
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associated  with  published  life  history  studies  are  summa* 
rized  in  Table  VII. 


TABLE  VII 

Suamary  of  Life  History  Threshold  Values 


|  ft 

« 

<5 


TIME  OF 

THRESHOLD 

VALUE 

STOTT 

LOCATION 

CASES 

TEAR 

INFRARED 

VISUAL 

Griffith, 

Florida, 

34  daps 

T1 - rf 

-20°C 

80  counts* 

•t  il„ 

Venezuela, 

1969-1976 

(1978) 

Honduras, 

and  hurricanes 

lapacting  Eaat 

Coast  United 

Stataa 

Stove, 

tropical  North 

37  obser¬ 

Saptaaber 

-26°C 

0.45  albedo 

•t  el.. 

Atlantic 

vations 

1974 

(sun  over¬ 

(1979) 

head) 

Wylie 

Montreal 

6  day* 

June  1977 

-16°C 

_ 

(1979) 

Saptaaber 

1977 

Negri  and 

OklihoMf 

1  day 

April  24, 

-27°C 

Adler 

Arkansas, 

(13  thunder-  1975 

(1981) 

Missouri 

atona) 

*  ATS-3  satellite 

\ 

Vk* 


The  Scofield/Oliver  (Scofield,  1981)  analysis  follows  a 
decision  tree  procedure  to  estimate  half-hourly  rainfall  for 
deep  convective  systeas  within  tropical  air  masses.  Osin? 
enhanced  infrared  and  high  resolution  visual  satellite 
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images,  the  technique  involves  first  identifying  the  active 
convective  portion  of  the  cloud,  or  cluster,  from  two  con¬ 
secutive  satellite  isages.  Once  the  active  portion  is  iden¬ 
tified,  the  half-hourly  rainfall  estimation  is  computed 
based  on  such  factors  as  cloud  top  temperature,  cloud 
growth,  and  departure  of  precipitable  water  from  a  summer¬ 
time  normal. 

The  Grif  fith/8oodley  {Griffith  et  al. ,  1978)  technique 

is  designed  to  estimate  rainfall  in  the  tropics,  over  large 
space  and  time  scales,  using  geosynchronous  visual  or  infra¬ 
red  satellite  imagery.  This  time-dependent  technique  was 
empirically  derived  as  a  relationship  between  cloud  area, 
echo  area,  and  rain  rate  for  two  areas  in  south  Florida, 
with  raingage-radar  providing  the  ground  truth,  and  was  then 
tested  in  other  tropical  areas.  This  scheme  was  subse¬ 
quently  tested  further  in  extra tropical  areas  (Griffith  et 
al.,  1980),  with  modifications  to  the  rainfall  amount 
predicted. 

The  determination  of  a  cloud  area-rainfall  relationship 
first  required  the  specification  of  both  a  visual  and  an 
infrared  threshold  to  define  the  cloud  area.  The  visual 
brightness  threshold,  normalized  for  radiation  geometry,  was 


80  counts  for  the  third  Application  Technology  Satellite 
(ATS -3)  and  the  infrared  threshold  was  253K  (-20°C)  .  The 
thresholds  were  based  on  a  comparison  of  the  clouds  with  a 
given  maximum  digital  count  and  the  radar  echoes  associated 
with  these  clouds. 

The  empirical  cloud  area -rainfall  relationship  was 
derived  as  a  two  step  process.  First,  a  relationship 
between  the  cloud  area  and  the  radar  echo  area,  normalized 
for  the  maximum  area  achieved  by  the  cloud  or  cluster,  was 
established  for  the  visible  and  infrared  satellite  data. 
Second,  the  relationship  between  the  echo  area  and  rain  vol¬ 
ume  was  determined  and  was  of  the  fora: 

Bv  =  I  Ae  (5) 

where  Rv  is  rain  volume  per  hour  (m3-h“M  ,  I  is  rain  in  units 
of  (m*-'km-*-h~t)  ,  and  Ae  is  the  echo  area  (km2)  defined  by 
the  1  mm-h”1  rain  rate.  Thus,  given  a  time  sequence  of  con¬ 
vective  clcuds  (or  cluster  areas)  measured  from  visible  or 
infrared  satellite  images,  volumetric  rain  rate  can  be 
estimated. 

Stout  ft  al.  (19  79)  modified  the  Griffith/Woodley  tech¬ 
nique  (Griffith  et  fl.,  1978)  to  estimate  volumetric  rain 
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rate  directly  from  a  cumulonimbus  cloud  area  and  area  change 
according  to  the  equation: 

B  =  aQA  ♦  at  dA/dt  (6) 

where  R  is  the  voluaetric  rainfall  of  the  cloud  (m3-s~l)»  A 
is  the  cloud  area  (a2)  ,  dA/dt  is  the  change  of  cloud  area 
over  tiae  (mz-s-1)#  and  a0  and  a,  are  constants  with  dimen¬ 
sions  m-s-*  and  a,  respectively.  The  two  constants  were 
calculated  by  a  least  squares  fit  of  cloud  area-rain  rate 
pairs  based  on  visible  and  infrared  geosynchronous  satellite 
data  and  5.3  ca  ship  radar  rain  data  collected  during  GATE. 
The  cloud  area  and  its  change  are  defined  by  the  threshold 
value.  The  visible  threshold  for  cloud  area  calculations 
was  60  digital  counts  on  the  ATS-3  (corresponding  to  an 
albedo  of  0.45  with  the  sun  overhead)  ,  or  equivalently  172 
digital  counts  on  the  first  Geosynchronous  Meteorological 
Satellite  (SMS  1).  The  infrared  threshold  was  160  digital 
counts  (-26°C).  The  standard  error  between  rhe  estimated 
rainfall  and  the  mean  radar  rainfall  was  62%  and  76%  for  the 
visual  and  infrared  equations  respectively. 

Wylie  (1979)  attempted  to  use  the  tropical  convective 
rainfall  techniques  of  Griffith  et  al. ,  (1978)  and  stout  et 


al . .  (1979)  for  estimating  precipitation  in  Montreal,  Can¬ 
ada.  Using  visual  satellite  data,  corrected  for  the  chang¬ 
ing  sun  angle  (Mosher,  1975),  infrared  satellite  data,  and 
10.0  cm  radar  meaured  rainfall  rates,  Wylie  studied  six  days 
of  precipitation,  three  days  each  in  June  and  September 
1977.  Wylie  concluded  that  because  of  air  mass  differences 
between  Montreal  and  the  tropics,  the  Griffith  and  Stout 
estimation  techniques  did  poorly  in  Montreal,  Canada.  The 
singlemost  important  limitation  with  these  two  schemes  was 
the  difficulty  of  measuring  cumulonimbus  cloud  area  when  the 
"anvils  were  often  merged  into  large  cloud  masses  and  the 
extensive  stratus  cloud  cover  often  obscured  the  pictures." 
Wylie  also  noted  that  the  Griffith  el.  al.  (1978)  threshold 
of  - 26°C  had  to  be  changed  to  -16°C  for  the  summertime  Mont¬ 
real,  Canada,  area.  With  the  warmer  cloud  top  temperatures 
the  cloud  areas  were  a  larger,  more  appropriate  size  for 
tracking. 

Wylie  (1979)  then  attempted  to  combine  sounding  data 
input  into  a  one- dimensional  model  (Simpson  and  wiggert, 
1969)  and  satellite  cloud  cover  measurements  to  estimate 
rainfall  for  Montreal.  With  the  GATE  measurements  for  rain 
rates  associated  with  sat ellite- derived  cloud  areas  and  the 


model  output,  rainfall  rates  were  estimated  by  multiplying 
the  two  values.  The  most  accurate  estimations  were  for  the 
cumulus  clouds  in  the  warm  air  masses  occurring  in  June,  the 
cases  the  model  was  designed  to  handle.  Wylie  concluded 
that  in  order  to  estimate  rainfall  in  all  geographical  areas 
and  seasons  a  more  sophisticated  model  would  be  needed. 

Negri  and  Adler  (1981)  did  one  case  study  of  fifteen 
thunderstorms  in  the  Oklahoma,  Arkansas,  and  Missouri  area 
on  24  April  1975.  They  used  radar  data  for  ground  truth  and 
had  special  5  minute  SOES-E  satellite  passes  over  the  area 
of  interest.  They  were  able  to  determine  that  the  precipi¬ 
tation  began  falling,  as  indicated  by  radar  data,  for  cloud 
top  temperatures  ranging  from  229K  to  260K  (-44°C  to  -13°C). 
The  mean  cloud  tcp  temperature  value  was  247K  (-26°C) . 
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III.  DATA  PROCESSING 

A.  INTRODUCTION 

The  data  set  asseabled  for  this  study  consists  of  collo¬ 
cated  GOES-E  satellite  data  and  Service-A  hourly  surface 
observations  for  the  southeastern  United  States  during 
August  1979.  The  GOES-E  data  consists  of  10  x  10  pixel 
matrices  of  visual  and  infrared  satellite  data  centered  over 
each  of  137  surface  stations  (Fig.  8)  all  south  of  40°N. 
The  satellite  data  are  measured  with  the  visual  Infrared 
Spin  Scanned  Radiometer  (VT.SSB)  which  have  subsatellite 
point  spatial  resolutions  of  1  and  7  km  for  the  visual  and 
infrared  channels,  respectively.  The  GOES-E  navigation  was 
completed  by  Ban-computer  Interactive  Data  Access  System 
(BcIDAS)  at  the  University  of  Wisconsin  using  the  full  reso¬ 
lution  visual  data,  with  an  accuracy  of  1-2  pixels  (1-2  km). 
The  full  resolution  visual  data  were  averaged  to  a  7  km  res¬ 
olution,  to  equal  the  infrared  data  resolution.  The  visual 


and  infrared  digital  counts  range  from  values  of  0-255.  The 
10  x  10  pixel  GOES-E  visual  and  infrared  satellite  data  each 
cover  an  area  45  nmi  x  45  nmi  at  30°N  (60  nmi  x  60  nmi  at 


Figure  8.  Geographical  Locations  of  Service-A  Station 
Report  Data 


The  data  are  divided  into  two  no-pracipitation  catego¬ 
ries  (Table  VIII)  and  seven  precipitation  categories  (Table 


k  . 


IX)  to  investigate  precipitation  specification,  convective 
versus  continuous  precipitation  specification,  and  qualita¬ 
tive  specification  of  light  versus  moderate/heavy  precipita¬ 
tion.  The  pixel  array  size  is  also  varied  from  the  10  x  10 
array  size  to  an  8  x  8,  a  6  x  6,  and  a  4  x  4  array  size  to 
investigate  the  differences  in  the  data  resulting  from  vari¬ 
ous  resolution  sizes  within  a  particular  weather  condition 
cl assification. 


B.  EAT A  SOBT 

For  the  combined  visual  and  infrared  threshold  specifi¬ 
cation  of  precipitation,  satellite  data  for  1200-2000  GHT, 
corresponding  to  0800-1600  EOT,  were  sorted  into  precipita¬ 
tion  and  no-precipitation  groups  (Fig.  9).  The  0800-1600 
EOT  interval  was  chosen  to  avoid  distortion  of  the  visual 
satellite  data  due  to  a  low  solar  elevation  angle.  The  vis¬ 
ual  data  were  normalized  and  converted  to  albedos  based  on 
the  work  done  by  Huench  and  Keegan  (1979).  This  scheme  cor¬ 
rects  for  the  varying  zenith  angle  as  well  as  adjusting  the 
visual  satellite  data  for  anisotropic  scattering  as  related 
to  the  zenith  angle.  (See  Appendix  A  for  further  specific 
information  concerning  the  Huench  and  Keegan  normalization 
scheme.) 


SURFACE  DATA 

East  Coast  and  Gulf  Coast  United 
States/  August  1979/  Service  -  A 
hourly  observations  <70/  623) 


COLLOCATED  SATELLITE  DATA 

Visual  10  x  10 

PIXEL  SIZE 

Infrared  10  x  10 

PIXEL  SIZE 

DAYLIGHT  TIME  CHECK 

(avoid  low  solar  angle  for 

VISUAL  SATELLITE  DATA) 

Include  1200-2000  GMT 
<0800-1600  EDT)  REPORTS 


Isolate  Applicable  Weather 
and  Cloud  Cases  for  Study 
(see  Tables  VI II  and  IX) 


Delete  reports  with  any 

ZERO  VALUES  IN  SATELLITE 
VISUAL  AND  INFRARED  DATA 


Normalize  Visual  Data 

CONVERT  VISUAL  DIGITAL 
COUNTS  TO  ALBEDOS 
(MUENCH  AND  KEEGAN/  1979) 


Compute  means  and  standard 

DEVIATIONS  OF  THE  ALBEDOS 
AND  CLOUD  TOP  TEMPERATURES 
FOR  EACH  10  X  10/  8  X  8,  6  X  6/ 
AND  4  X  4  PIXEL  ARRAY  SIZE 


CALCULATE  MEAN  AND  STANDARD 
DEVIATION  OF  THE  MEANS  AND 
STANDARD  DEVIATIONS  FOR  EACH 
WEATHER/CLOUD  CLASSIFICATION 
AND  PIXEL  ARRAY  SIZE 


Figure  9.  Flow  Chart  of  Data  Processing 
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While  albedo  values  cannot  exceed  1.00,  the  Muench  and 


Keegan  (1979)  scheme  allows  the  values  to  overshoot  1.00,  up 
to  a  value  of  1.20.  Therefore,  the  visual  satellite  values 
are  not  true  albedos,  but  estimated  albedos.  The  extended 
visual  noraalized  data  scale  was  used  to  facilitate  compari¬ 
son  of  the  results  in  this  effort  to  the  most  extensive  bi- 
spectral  threshold  precipitation  specification  of  Muench  and 
Keegan  (1979).  The  Muench  and  Keegan  (1979)  normalization 
scheme  specifies  that  any  computed  albedo  greater  than  1.20 
be  set  equal  to  1.20  to  limit  the  unreasonably  large  values. 
Similiarly,  the  scheme  specifies  computed  albedos  less  than 
0.15  be  interpreted  as  the  ground  or  water  surface  reflec¬ 
tance  and  the  value  0.00  be  assigned.  The  infrared  data 
were  processed  in  digital  counts  and  converted  to  cloud  top 
temperatures  prior  to  statistical  computations  and  graphical 
displays. 

The  no-precipitation  data  (Table  VIII)  are  comprised  of 
the  digital  visual  and  infrared  10  x  10  pixel  arrays  of 
those  Service-i  station  reports  not  showing  any  "RM  in  the 
current  weather  group.  Thus,  the  no-precipitation  group 
includes  stations  reporting  drizzle  (weather  codes  L-,  L, 
and  L4) . 
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The  no-precipitation  data  are  divided  into  two  catego¬ 
ries,  overcast  ceiling  (category  2A)  and  overcast  and  broken 
ceiling  (category  2B)  .  Cloud  cover  is  based  on  the  three 
digit  cloud  group  in  the  Service -A  surface  observation.  The 
first  digit  indicates  the  amount  of  low  clouds,  where  0  is 
defined  as  clear,  1  is  scattered  (one-eighth  to  four- 
eighths  cloud  cover),  2  is  broken  (five-eighths  to  seven- 
eighths.  cloud  cover),  and  3  is  overcast  (eight-eighths  cloud 
cover) .  The  second  and  third  digit  indicate  the  amount  of 
■iddle  and  high  clouds,  respectively.  The  sane  0-3  values 
defined  for  the  low  clouds  are  used  for  middle  and  high 
cloud  aaount. 
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The  precipitation  data  (Table  IX)  are  comprised  of  the 


visual  and  infrared  10  x  10  pixel  arrays  of  those  Service-A 


station  reports  showing  any  *•  RM  in  the  current  weather 


group.  Two  precipitation  observations  were  excluded  from 


the  data  set  because  each  report  also  indicated  clear  skies. 


TABLE  IX 

Classification  of  Precipitation  Data  Groups 


Service- A 

Category  Same  Current  Weather  Reports 


1 

Precipitation 

any  "R" 

5  38 

1  A 

Precipitation 
Overcast  Ceiling 

any  "R"  and  overcast 
ceiling  (as  defined  in 
Table  vlll  category  2A) 

329 

IB 

Precipitation 
Overcast  and 
Broken  ceiling 

any  "R"  and  overcast 
ana  broken  ceiling  (as 
defined  in  Table  VIII 
category  2B) 

534 

1C 

Continuous 

Precipitation 

R-  ,  B,  R+ 

1 12 

ID 

Convective 

Precipitation 

RW-,  R W,  RW+,TRW-, 

TRW,  TRW+,  TR-,  TR,  TR 

426 

IE 

Light 

Precipitation 

R-,  RW-,  TRW-,  TR- 

464 

IP 

Moderate/Heavy 

Precipitation 

R,  R  +  ,  RW,  RW  +  ,  TRW, 
TRW+ ,  TR,  TR 

74 

The  general  precipitation  data  (category  1)  are  divided 
into  six  groups:  precipitation  overcast  ceiling  (category 


I A )  ,  precipitation  overcast  and  broken  ceiling  (category 

IB) ,  continuous  (category  1C),  convective  (category  ID), 
light  (category  IE),  and  aodera  te/heavy  (category  IF)  pre¬ 
cipitation.  These  seven  precipitation  groups  are  used  to 
investigate  precipitation  specification,  convective  versus 
continuous  precipitation  specification,  and  qualitative 
specification  of  light  versus  mo  derat  e/heavy  precipitation. 

C.  STATISTICAL  TREATMENT 

The  means  and  standard  deviations  of  albedos  and  cloud 
top  temperatures  of  each  10  z  10  pixel  array  for  the  weather 
types  listed  in  Tables  VIII  and  IX  were  calculated.  Means 
and  standard  deviations  of  albedos  and  cloud  top  tempera¬ 
tures  were  also  calculated  for  the  8x8,  6x6,  and  4x4 

pixel  arrays  centered  over  the  surface  station.  The  8x8, 
6x6,  and  4x4  pixel  arrays  are  equal  to  36  nmi  x  36  nmi, 
27  nmi  x  27  nmi,  and  22  nmi  x  22  nmi  at  30°N  respectively. 
Variation  of  the  digital  satellite  areal  coverage  is  used  to 
investigate  the  differences  in  the  statistics  due  to  the 
chosen  resolution  size. 

The  data  sets  in  Tables  VIII  and  IX  are  represented, 
first,  by  the  mean  and  standard  deviation  of  the  resolution 
cell  means  and  standard  deviations.  Second,  these  data  sets 
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are  represented  by  the  distributions  of  the  mean  cloud  top 
temperatures  and  albedos  where  the  nean  cloud  top  tempera¬ 
tures  are  sorted  into  ten  Kelvin  (K)  intervals  and  the  mean 
albedos  are  sorted  into  0.10  intervals.  These  representa¬ 
tive  statistics  and  distributions  are  calculated  for  the 


four  pixel  array  sizes. 

The  statistical  and  distribution  results  for  differing 
resolution  sizes#  bi-spectral  threshold  specification  of 
precipitation#  and  separation  of  light  fron  moderate/heavy 
precipitation  are  discussed  in  Chapter  IT. 
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IV.  8ES0LTS 


A.  INTRO D OCT ION 

The  figures  presented  in  this  chapter  display  the  dis¬ 
tributions  of  the  grand  leans  of  the  resolution  cell  leans 
of  albedos  and  cloud  top  temperatures  for  the  data  sets 
listed  in  Tables  VIII  and  IX  for  the  four  array  sizes.  The 
■ean  values  are  sorted  into  ten  Kelvin  intervals  and  0.10 
estimated  albedo  intervals. 

B.  RESOLUTION 

The  effect  of  satellite  resolution  in  representing  gen¬ 
eral  precipitation  (category  1)  and  no-precipitation  over¬ 
cast  (category  2A)  data  are  explored  for  four  resolution 
sizes.  The  four  sizes  are  10  x  10,  8  x  8,  6  x  6,  and  4  x  u 
and  are  approximately  equal  to  areas  of  2025  nmi2,  1296 
nmi2,  729  nmi2,  and  484  n«i2  at  30°N  respectively. 

The  general  precipitation  (category  1)  and  no-precipita¬ 
tion  overcast  (category  2A)  data  were  chosen  for  study 
because,  while  they  represent  two  different  weather  condi¬ 
tions,  their  albedo  and  cloud  top  temperature  distributions 
have  the  largest  amount  cf  overlap  when  compared  to  any 


51 


VvvV'-.-'v'l 


other  pair  of  precipitation  versus  no-precipitation  data 
sets.  The  possibility  arises  that  statistical  differences 
in  the  four  resolution  sizes  might  be  sufficient  or  comple¬ 
ment  other  information  in  delineating  these  two  weather 
conditions. 

1.  Precipitation  Data 

The  general  precipitation  (category  1)  data  are  com¬ 
prised  of  329  overcast  ceiling  reports  (61%),  205  broken 

ceiling  reports  (38%),  and  4  scattered  ceiling  reports  (1%). 

a.  Bean  Statistics 

The  pracipitation  data  (category  1)  differences 
between  the  means  of  the  cell  means  visual  and  infrared  10  x 
10  and  4  x  4  array  sizes  are  0.  035  and  2.2K,  respectively 
(Table  X)  .  The  trend  of  the  mean  of  the  means  is  toward 
higher  albedo  values  and  colder  cloud  top  temperatures  with 
the  decreasing  area  or  array  size.  The  standard  deviations 
of  the  means  siocLliarly  show  an  increase  in  the  albedo, 
0.016,  and  cloud  top  temperature,  1.0K,  from  the  10  x  10 
array  size  to  the  4x4  array  size. 

b.  Standard  Deviation  Statistics 

The  standard  deviation  statistics  display  the 
opposite  trend  with  decreasing  area  as  the  mean  statistics. 


c.  Distribution  Discussion 

The  distributions  of  the  precipitation  data  are 
shown  in  Pigs.  10,  11,  12,  and  13.  There  is  a  discernible 
shift  toward  higher  albedcs  and  colder  cloud  top  tempera¬ 
tures  of  the  2%  and  3%  frequency  isopleth  with  decreasing 
array  size.  This  upward  shift  is  also  reflected  in  the  mean 
of  the  means  (Table  X).  The  appearance  of  the  5%  frequency 
isopleth  in  the  6x6  and  4  x  4  array  sizes  at  high  albedos 
and  cold  cloud  tcp  temperatures  highlights  the  shift. 
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Figure  10. 
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Precipitation  Data  for  10  x  10  Array  Size  (The 
mean,  .579  and  253. OK,  interval  is  boxed.  The 
2%  and  3*  frequencies  are  for  11  and  16 
occurrences,  respectively.) 
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Figure  13.  Precipitation  Data  for  4x4  Array  size  (The 

mean.  .614  and  250. 8K,  interval  is  boxed.  The 
2S,  3%,  and  5%  frequencies  are  for  11,  16,  and 
27  occurrences,  respectively.) 

The  distr  ibutions  of  the  10  x  10  and  8x8  array 
sizes  (Figs.  10  and  11)  are  unimodal  while  the  6x6  and  4  x 
4  array  sizes  (Figs.  12  and  13)  appear  to  be  more  bimodal. 
The  four  array  sizes  were  tested  for  a  Gaussian  distribution 
with  the  Chi-square  test  and  all  failed  at  any  confidence 
level.  Therefore,  differences  in  the  four  resolutions  can¬ 
not  be  adequately  tested  by  well  defined  statistical  methods 
based  on  an  assumed  normal  distribution. 

The  similiarities  between  the  10  x  10  and  8x8 


array  sizes  (Figs.  10  and  11)  and  the  6x6  and  4x4  array 


sizes  (Pigs.  12  and  13)  are  further  illustrated  in  Pig.  14. 
k  diagonal  cut  is  plotted  for  each  of  the  four  array  sizes 
where  the  lines  plotted  are  shown  as  a  dashed  boxed  area  in 
Pigs.  10-13.  Tha  diagonal  cut  reveals  the  close  agreement 
between  the  6x6  and  4x4  array  sizes  along  the  line.  The 
10  x  10  and  8x8  array  size  lines  follow  the  same  general 
trend  but  do  not  coincide  as  closely  as  tha  6x6  and  4x4 
array  size  lines. 

The  chosen  diagonal  line  results  in  the  8  x  8 
array  size  distribution  appearing  more  smoothed  than  the  10 
x  10  (Pig.  14),  as  there  is  no  relative  minima  at  the  4.0 
interval  for  the  8x8  array  size.  Tha  10  x  10  array  size, 
with  the  greater  areal  extent  and  therefore  more  averaging 
of  differing  clouds  and  clear  areas,  is  expected  to  possess 
the  "smoo thest"  appearance,  the  lowest  number  of  relative 
maxima  and  minima  of  the  four  array  sizes.  However,  the  8  x 
8  array  size  actually  displays  the  fewest  relative  maxima 
and  minima  along  the  chosen  diagonal  line.  The  smoother  8  x 
8  array  size  cannot  be  explained  in  terms  of  significant 
differences  in  the  number  of  cases  of  different  ceiling 
types  or  different  precipitation  types  occurring  in  interval 
4.0  between  the  four  array  sizes.  Quite  simply,  the 


smoother  8x8  array  size  apparently  results  from  the  sort 
ing  intervals  chosen  for  the  distributions. 


1.0  2.0  3.0  4.0  5.0  6.0  7.0  8.0  9.0  10.0  11.0 


CLOUD  TOP  TEMPERATURE/  ALBEDO  INTERVALS 


Figure  14.  Precipitation  Array  Size  Distributions  Along 

Diagonal  Line  (The  line  represents  the  10  x  10, 
dotted  line  the  8  x  8,  dashed  line  the  6x6, 
and  dash  dot  line  the  4x4  array  size.) 


The  elongated  shape  of  all  four  precipitation 


distribution  array  sizes  (Figs.  10,  11,  12,  and  13)  reveal 


the  variation  in  the  areal  amount  of  cloudiness  and  precipi¬ 


tation.  The  distributions  range  from  high  albedos  and  cold 


cloud  top  temperatures  (indicative  of  satellite  fields  of 


view  filled  with  precipitating  clouds)  to  low  albedos  and 
warn  cloud  top  temperatures  (indicative  of  satellite  fields 
of  view  partially  filled  with  precipitating  clouds).  The 
visual  and  infrared  satellite  data  distrib unions  in  Figs. 
10 ,  11,  12,  and  13  agree  with  the  elongated  shapes  of  Platt 
(1981)  for  his  cloud  cl assifications  of  cumulus,  frontal, 
and  jetstream  cirrus  clouds  and  agree  with  Coakley  and 
Bretherton  (1982)  for  their  general  clouds  present  in  a  1000 
Jem2  Pacific  Ocean  area, 
d.  Summary 

&  satellite  field  of  view  filled  with  a  precipi¬ 
tating  cloud  is  expected  to  have  higher  albedo  and  colder 
cloud  top  temperature  values  than  a  partially  filled  field 
of  view.  Additionally,  the  filled  field  of  view  would  have 
a  more  uniform  texture,  as  reflected  in  variance  or  standard 
deviation  values,  than  a  partially  filled  field  of  view. 
The  statistics  discussed  in  this  study  confirm  these  expec¬ 
tations  for  this  data  set.  As  the  array  size  or  field  of 
view  is  decreased,  the  mean  statistics  increase  while  the 
standard  deviation  statistics  decrease  (Table  X) . 

Fcr  the  precipitation  data  (category  1)  ,  there 


are  significant  differences  between  the  four  resolution 


sizes.  These  differences  are  reflected  in  the  upward  trend 
in  albedos  and  cclder  cloud  top  temperatures  of  the  mean  and 
the  standard  deviation  of  the  Beans  with  decreasing  area  or 
array  size  (Table  X)  .  The  reverse  trend  is  found  in  the 
mean  and  standard  deviation  of  the  standard  deviations. 

Tke  distributions  have  significant  differences 
also.  The  relatively  ccarse  resolution  10  z  10  and  8x8 
array  sizes  have  a  uniaodal  distribution  while  the  rela¬ 
tively  fine  resclution  6  z  6  and  4  z  4  have  a  bimodal 
di stribution. 

The  four  array  sizes  discussed  vary  in  their 
statistics  and  distributions  in  representing  the  precipita¬ 
tion  (category  1)  data.  The  choice  of  satellite  resolution 
for  representation  of  the  precipitation  data  will  influence 
comparison  of  these  data  with  other  data.  Therefore,  for 
the  remainder  of  this  study  the  precipitation  data  for  all 
classifications  will  be  discussed  using  both  the  10  x  10  and 
4x4  array  size.  , 

One  additional  topic  to  explore  is  that  a  number 
of  cases  with  albedos  less  than  0.40  appear  in  the  precipi¬ 
tation  data  (category  1)  in  all  four  size  distributions 
(Figs.  10,  11,  12,  and  13).  These  low  albedo  values  suggest 
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the  possibility  that  there  eight  be  a  consistent  low  bias  in 
the  normalization  scheme.  But  a  closer  look  at  the 
individual  reports  with  albedos  less  than  0.40  show  no 
pattern  involved  with  either  the  GHT  hour  or  the  longitude 
or  latitude  of  these  station  reports.  Further  analysis  of 
these  low  albedo  precipitation  reports  are  discussed  in  the 
light  precipitation  section  (IV. E.)  . 

2.  Hp~pregipitation  Overcast  Data 
a.  Mean  Statistics 

For  the  no-precipitation  overcast  cases  (cat¬ 
egory  2k)  ,  differences  between  the  means  of  the  means  visual 
and  infrared  10  x  10  and  4x4  array  sizes  are  0.011  and  0.3 
K,  respectively  (Table  XI).  The  0.3K  infrared  difference  is 
within  the  0.5K  noise  level  of  the  VISSR  infrared  sensor. 
The  trend  of  the  visual  mean  of  the  means  is  upward  with  the 
decreasing  array  size.  The  standard  deviations  of  the  means 
show  an  increase  in  the  albedo,  0.014,  and  cloud  top  temper¬ 
ature,  0.8K,  from  the  10  x  10  to  4  x  4  array  sizes.  Once 
again,  both  mean  statistics  have  an  upward  trend  with 
decreasing  array  size. 
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TABLE  XI 

Mo-Precipitation  Data  Statistics  for  Pour  Array  Sizes 

Overcast  Ceiling 


11  I  11 

ill 

6x6 

4x4 

Mean  of  (VIS) 

Means 

.  407 

.411 

.415 

.418 

<IR) 

272.6K 

272. 5K 

272. 4K 

272. 3K 

(-13C) 

(-  1°C) 

<-1°C) 

(-1°C) 

Standard  (VIS) 
Deviations 

.  205 

.  208 

.213 

.219 

of  Means  (IB) 

1  7.  8K 

1  8.  IK 

18.3K 

18. 6K 

Mean  of  (VIS) 

Standard 

.  122 

.  113 

.102 

.  087 

De  viations  (IB) 

Standard 

5.4K 

4.7K 

3 . 9K 

2.9K 

Deviation  (VIS) 
of  the 

.056 

.  054 

.052 

.049 

Standard 
Deviations  (IB) 

4.  8K 

4.4K 

3 . 9K 

3.  IK 

b.  Standard  Deviation  Statistics 

Conversely,  the 

standard 

deviation 

statistics 

have  a  downward 

trend  with 

deer  easing 

area  size. 

The  dif- 

ferences  between 

the  10  x 

10  and  4 

x  4  array 

visual  and 

infrared  aeans  of  the  standard  deviations  are  0.035  and 


2.5K,  respectively.  The  differences  in  the  standard  devia¬ 
tions  of  the  standard  deviations  are  0.007  in  the  visual  and 
1.7K  in  the  infrared  values. 


c.  Distribution  Discussion 

Tbs  distributions  of  the  no-precipitation  over¬ 
cast  data  are  shore  in  Pigs.  15,  16,  17,  and  18.  These  dis¬ 
tributions  are  quite  different  from  the  precipitation 
distributions  (Figs.  10,  11,  12,  and  13).  As  expected,  the 
no-precipitaticn  overcast  data  are  clustered  at  the  low 
albedo  and  warm  cloud  top  temperature  values.  The  5%  fre¬ 
quency  isopleths  in  Figs.  15,  16,  and  17  show  a  grouping  of 
the  data  at  albedos  ranging  from  0.30  to  0.50  and  cloud  top 
temperatures  from  280K  to  290K.  A  bimodal  distribution 
appears  in  the  finer  resolution  8x8,  6x6,  and  4x4 

array  sizes  (Figs.  16,  17,  and  18).  This  shifting  of  the 
no-precipitaticn  overcast  data  into  two  relative  maxima  for 
the  three  smallest  array  sizes  is  the  sole  significant  dif¬ 
ference  in  the  four  distributions. 

d.  Summary 

The  no-precipitation  overcast  data  display  an 
upward  trend  in  the  two  mean  statistics  with  decreasing 
array  size,  although  the  0. 3K  mean  of  the  means  infrared 
difference  is  not  significant.  Conversely  the  two  standard 
deviation  statistics  decrease  with  decreasing  array  size. 
These  trends  are  consistent  with  the  expected  statistical 


CLOUD  TOP  TEMPERATURE  (X) 

Figure  17.  No-Precipitation  Overcast  Data  for  6x6  array 
Site  (The  lean,  .415  and  272.  4K,  interval  is 
boxed.  The  256  ,  3  56,  and  5%  frequencies  are  for 
40,  59,  and  99  occurrences,  respectively.) 


CLOUD  TOP  TEMPEXATUIK  (X)  , 

Figure  18.  No-Precipitation  Overcast  Data  for  4x4  array 
Size  (The  lean,  .418  and  272. 3K,  interval  is 
boxed.  The  2 X ,  356,  and  5  56  frequencies  are  for 
40,  59,  and  99  occurrences,  respectively.) 


.'.V, 


trends  discussed  in  the  precipitation  section  (IV.B.1.). 
However,  in  these  data,  there  is  a  greater  similiarity  in 
the  statistics  for  the  four  sizes  because  the  ceilings  are 
all  overcast  reports. 

The  distributions  (Figs.  15,  16,  17,  and  18) 
provide  visual  confirmation  of  the  similarities  between 
each  of  the  four  array  sizes.  With  the  exception  of  the 
second  relative  maxima  at  albedos  of  0.30  to  0.40  and  cloud 
top  temperatures  of  260K  tc  270K  appearing  in  the  8x8,  6  x 
6,  and  4x4  array  sizes  (Figs.  16,  17,  and  18),  the  four 
distributions  are  nearly  identical.  Because  there  are  sta¬ 
tistical  differences  and  distributional  differences  in  the 
four  array  sizes,  the  10  x  10  and  4x4  arrays  sizes  will  be 
used  to  represent  the  two  no-precipitation  data  categories. 

3.  Precipitation  and  Ho- precipitation  comparison 

The  possibility  arises  that  statistical  differences 
in  the  four  resolution  sizes  might  be  sufficient,  or  at 
least  complement  other  information,  in  delineating  the  pre¬ 
cipitation  frcm  the  nc- precipita tion  weather  condition.  The 
guestion  then  becomes,  is  there  a  statistic  associated  with 
variation  of  the  array  size  within  Tables  X  and  XI  which 
differentiates  the  precipitation  reports  from  the  no-precip- 
itation  overcast  reports? 
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If  the  satellite  data  processor  can  vary  the 
resolution  size,  as  was  dene  in  this  study  by  simply  averag¬ 
ing  different  pixel  array  sizes,  a  trend  in  the  visual  and 
infrared  data  might  be  used  to  differentiate  rhese  two 
weather  conditions.  The  most  significant  trend  difference 
in  the  precipitation  data  (Table  X)  and  no-precipitation 
overcast  data  (Table  XI)  occurs  in  the  mean  of  the  means. 
Recall  that  the  mean  of  the  means  precipitation  visual  dif¬ 
ference  between  the  10  x  10  and  4x4  array  sizes  was  0.035 
while  the  no -precipitation  overcast  visual  difference  was 
0.011.  Similiarly  the  infrared  differences  were  2.2K  and 
0.3K  for  the  precipitation  and  no-precipitation  overcast, 
respectively.  The  precipitation  data  show  a  greater  upward 
trend  toward  higher  albedos  and  colder  cloud  top  tempera¬ 
tures  than  the  no-precipitation  overcast  with  the  finer  sat¬ 
ellite  resolution. 

Differentiation  between  these  two  data  sets  based  on 
a  comparison  of  the  trend  in  the  mean  of  the  means  are  sug¬ 
gested  by  Tables  X  and  XI.  It  must  be  emphasized  that  these 
tables  are  based  on  many  reports  and  therefore  reflect  the 
most  typical  values.  Individual  reports  within  a  given 
interval  should  be  studied  to  provide  conclusive  evidence  as 
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to  whether  these  statistics  can  be  used  on  a  few  reports  to 
differentiate  precipitation  from  no-precipitation  overcast 
re  ports. 

C.  PBECI PIT ATICN  SPECIFICATION 

An  essential  difference  between  this  specification  study 
and  most  of  those  in  the  literature  is  that  the  distribu¬ 
tions  of  the  precipitation  and  no-precipitation  data  sets 
are  examined  in  detail  to  extract  information  about  the 
probability  of  correct  classifications.  Only  Lovejoy  and 
Austin  (1979)  present  their  data  distributions.  The  bi- 
spectral  and  life  history  method  thresholds  (Tables  I  and 
VII)  refer  to  the  typical  or  most  common  threshold  values 
for  precipitation,  which  is  assumed  to  be  equivalent  to  the 
mean  of  the  means  in  this  study.  Therefore,  the  mean  of  the 
means  can  be  compared  to  the  threshold  values  in  Tables  I 
and  VII.  Additionally,  a  bi -spectral  threshold  can  be  pro¬ 
posed  based  on  the  distributions  and  with  these  distribu¬ 
tions  the  amount  of  overlap,  or  the  percentage  of  correctly 
classified  precipitation  or  no-precipitation  cases,  can  be 


calculated . 


1  •  Overcast  Ceilings 

Are  the  precipitation  overcast  (category  1A)  and 
no-precipitation  overcast  (category  2A)  data  sets  suffi¬ 
ciently  separated  to  allow  differentiation  of  the  two  popu¬ 
lations?  If  so,  how  ouch  overlap  is  there  between  the  two 
data  sets? 

a.  Mean  of  the  deans 

The  mean  of  the  means  statistics  for  the  precip¬ 
itation  overcast  versus  no-precipitation  overcast  data. 
Table  XII,  shew  there  is  a  .242  and  24. 4K  difference  between 
the  two  populations  for  the  10  x  10  array  size  and  a  .254 
and  25. 9K  difference  for  the  4x4  array  size.  The  respec¬ 
tive  differences  are  greater  than  one  standard  deviation  of 
the  means  of  either  of  the  two  populations. 

If  the  two  precipitation  overcast  array  sizes 
are  compared  to  the  no- precipitation  overcast  array  sizes, 
the  stronger  trend  in  the  mean  of  the  means  (Table  XII)  is 
seen  in  the  precipitation  overcast  data.  While  the  precipi¬ 
tation  visual  and  infrared  values  vary  by  0.023  and  1.8K 
between  the  10  x  10  and  4x4  array  sizes,  the  no- precipita¬ 
tion  visual  and  infrared  values  vary  by  0.011  and  0.3K. 
Both  the  mean  cf  the  means  and  their  trends  can  be  used  to 


TABLE  XII 

Precipitation  Specification  Overcast  ceilings 


dean  of 

(VIS) 

Precipi¬ 

tation 

1.0  x  10 

.  64  9 

No -Preci¬ 
pitation 

10  x  10 

.  407 

Precipi¬ 

tation 

4x4 

.672 

No-Preci¬ 

pitation 

a  i  a 

.418 

deans 

(IR) 

24  8.2K 

272.  6K 

246. 4K 

272. 3K 

Standard 

(VIS) 

(-  25°C) 

.  186 

(-  1°C) 

.  205 

(-27°C) 

.201 

<-1°C) 

.219 

Deviations 
of  deans 

(IR) 

18. 9K 

1  7.8K 

19. 8K 

18. 6K 

dean  of 

(VIS) 

.  152 

.  122 

.108 

.  087 

Standard 

De  viations  (IB) 

8.  9  K 

5.4K 

4 . 4K 

2.9K 

Standard 

De  viation 

(VIS) 

.078 

.  056 

.069 

.  049 

of  the 

Standard 
Deviations  (IR) 

6.3K 

4.8K 

4.  OK 

3.  IK 

differentiate  precipitation  from  no-precipitation  for  a 
large  number  of  reports  in  a  region  similiar  to  this  summer¬ 
time  convective  shower  dominated  area. 

b.  Mean  of  the  Standard  Deviations 

The  visual  and  infrared  means  of  the  standard 
deviations  vary  by  0.  030  and  3.5K  for  the  10  x  10  array  size 
and  by  0.021  and  1.5K  for  the  4  x  4  nrray  size.  In  the 
infrared  values,  the  no-precipitation  mean  of  the  standard 
deviations  have  a  magnitude  60%  of  the  precipitation  values 
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and  produce  relatively  large  differences.  In  the  visual 
values,  the  no-precipitaticn  nea  n  of  the  standard  deviations 
have  a  magnitude  80%  of  the  precipitation  values.  The  rela¬ 
tively  large  differences  (only  the  mean  of  the  means  have  a 
larger  difference)  suggest  the  use  of  this  statistic  to  dif¬ 
ferentiate  precipitation  overcast  from  no- precipitation 
overcast. 

c.  Standard  Deviation  of  the  Means 

The  differences  in  the  standard  deviations  of 
the  means  (Table  XII)  are  nearly  egual  when  comparing  the 
two  10  x  10  array  sizes  and  the  two  4x4  array  sizes.  The 
visual  differences  are  0.019  for  the  10  x  10  and  0.018  for 
the  4x4  array  size.  Similiarly,  the  infrared  differences 
are  1.1K  for  the  10  x  10  and  1.2K  for  the  4x4  array  size. 
The  differences  in  the  4x4  array  size  (0.018  and  1.2K)  are 
comparable  to  the  relatively  significant  differences  in  the 
means  of  the  standard  deviations  (.021  and  1.5K).  However, 
the  differences  are  not  comparable  in  the  10  x  10  array 
size . 

d.  Standard  Deviation  of  the  Standard  Deviations 
The  differences  in  the  standard  deviations  of 

the  standard  deviations  (Table  XII)  for  the  two  array  sizes 


are  approximately  equal  also.  The  visual  differences  are 
0.022  for  the  10  x  10  and  0.  020  for  the  4  x  4  array  size. 
The  infrared  differences  are  1.5K  for  the  10  x  10  and  0.9K 
for  the  4x4  array  size.  Once  again  the  differences  in  the 
4x4  array  size  (0.020  and  0.9K)  are  approximately  equal  to 
the  differences  in  the  means  of  the  standard  deviations 
(0.021  and  1.5K),  particularly  in  the  visual  value, 
e.  Distribution  Discussion 

The  distributions  for  the  10  x  10  array  size 
precipitation  overcast  (Fig.  19)  and  no-precipitation  over¬ 
cast  (Fig.  15)  and  the  4x4  array  size  precipitation  over¬ 
cast  (Fig.  20)  and  no-Frecipi tation  overcast  ((Fig.  18) 
allow  visual  confirmation  of  the  degree  of  separation  of 
these  two  populations.  These  figures  verify  the  separation 
between  the  the  occurrence  maxima  of  the  two  populations 
while  showing  that  there  is  overlap  of  some  of  the  values  in 
the  two  populations. 

The  appearance  of  a  bimodal  distribution  in  the 
relatively  fine  resolution  4  x  4  array  size  precipitation 
overcast  (category  1A)  and  no-precipitation  overcast  (cat¬ 
egory  2&)  in  Figs.  20  and  18  cannot  be  explained  in  terms  of 
the  precipitation  categories  (Table  IX)  or  ceiling  cover. 


An  alternate  explanation  night  be  that  these  two  maxima 
reflect  different  synoptic  signatures.  Pour  weak  frontal 
systems  impact  this  data  set  region  during  August  1979  and 
cause  a  change  in  the  cloudiness  and  precipitation  pattern 
which  is  normally  produced  by  daytime  heating.  Recall  that 
the  Love  joy  and  Austin  (1979)  data  set  also  showed  bimodal 
distributions  for  the  cumulus  no-rain  reports.  The  bimodal 
distributions  in  this  study  may  not  be  due  to  a  synoptic 
signature.  Nonetheless,  this  possibility  should  be 
investigated. 

f.  Precipitation  Probabilities 

Precipitation  probabilities  (Figs.  21  and  22) 
were  computed  from  the  precipitation  overcast  (category  1A) 
and  the  no-precipitation  overcast  (category  2A)  data  for  the 
10  x  10  (Figs.  19  and  15)  and  the  4x4  (Figs.  20  and  18) 
array  sizes.  Estimated  albedos  greater  than  1.00  were  not 
included  in  these  probabilities  as  they  accounted  for  only 
two  and  three  no- precipitation  overcast  reports  and  four  and 
five  precipitation  overcast  reports  in  the  10  x  10  size  and 
4x4  array  size,  respectively.  The  two  probability  figures 
indicate  that  the  50%  probability  of  precipitation  is  not  a 
simple  function  of  mean  albedo  and  mean  cloud  top 
temperature. 
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The  50%  probability  line  (Fig.  21)  shows  the 
precipitation  at  low  clcud  top  temperatures  (270K-290K) 
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occurs  at  high  albedos  (0.80-1.  00)  and  at  cold  cloud  top 
temperatures  (210K-230K)  occurs  at  relatively  low  albedos 
(0.40-0.60).  One  exception  occurs  in  the  0.60-0.80  albedo 
and  210K-230K  cloud  top  temperature  interval  and  represents 
26  precipitation  reports  cf  64  total  reports.  The  100%  pre¬ 
cipitation  probability  at  0.00-0.20  albedo  and  230K-250K 
results  from  two  precipitation  reports.  Two  reports  in  a 
0.20  albedc  and  20K  cloud  top  temperature  interval  are  not  a 
sufficient  number  of  reports  to  be  a  significant  indication 
of  a  high  probability  of  precipitation.  The  50%  precipita¬ 
tion  probability  line  (Pig.  22)  in  the  4x4  array  size  data 
shows  the  same  general  trend  of  decreasing  albedo  with 
decreasing  cloud  top  temperatures.  In  this  fine  resolution 
data  (Pig.  22),  there  is  an  upturning  of  the  50%  probability 
line  at  the  coldest  cloud  top  temperatures,  210K-230K. 
There  are  73  no-precipitaticn  overcast  and  80  precipitation 
overcast  reports  in  the  4x4  array  size  data  in  the 
210K-230K  interval  so  the  upturning  is  not  the  result  of 
lack  of  data.  The  appearance,  once  again,  of  a  greater  than 
50%  probability  cf  precipitation  at  low  albedos,  0.00-0.20, 
between  210K-250K  results  from  a  total  of  four  reports  (one 
precipitation  report  of  two  total  reports  in  the  230K-250K 
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interval  and  two  precipitation  reports  of  two  reports  in  the 
210K-230K)  and  is  not  a  significant  indication  of  high  pre¬ 
cipitation  probability. 

If  a  straight  line  is  drawn  to  represent  the  50% 
probability  line,  the  lines  for  the  two  array  sizes  are 
nearly  coincident  until  the  230 K  cloud  top  temperature  is 
reached.  For  the  purposes  of  this  study,  a  simple  linear 
function  bi-spectral  precipitation  threshold  based  on  a  50% 
probability  of  precipitation  can  be  approximately  defined  as 
extending  from  1.00  albedo  and  2  90K  cloud  top  temperature  to 
0.60  albedo  and  210K  cloud  top  temperature. 

Comparison  of  this  linear  bi-spectral  threshold 
with  the  Huench  and  Keegan  (1979)  threshold  (Fig.  2),  shows 
the  proposed  threshold  has  an  albedo  approximately  0.10 
smaller  at  corresponding  cloud  top  temperatures.  The  50% 
probability  asymptote  at  the  warm  cloud  top  temperatures 
shown  in  the  Huench  and  Keegan  (  1979)  results  (Fig.  2),  are 
not  shown  in  Figs.  21  and  22  due  to  lack  of  reports  in  these 
values.  The  lower  albedo  values  for  this  linear  bi-spectral 
threshold  may  be  associated  with  the  dominance  of  convective 
precipitation  (426  of  5  38  reports)  in  the  precipitation 
data.  The  effect  of  convective  precipitation  reports  on  the 
satellite  albedo  values  is  discussed  in  Section  IV. D. 


g.  Suninary 


The  two  significant  statistics  for  precipitation 
specification  of  overcast  ceiling  reports  for  both  array 
sizes  are  the  mean  of  the  means  and  the  mean  of  the  standard 
deviations  (Table  XII).  The  probabilites  in  Figs.  21  and  22 
make  use  of  the  mean  of  the  means  only.  Inclusion  of  the 
mean  of  the  standard  deviations  in  the  precipitation  and 
ao-precipitaticn  distribution  plots  may  more  distinctly 
define  the  two  weather  data  types.  The  question  is,  how  to 
graphically  display  four  variables  (i.e.  four  dimensions)  in 
one  plot? 

One  solution  is  to  find  a  three-dimensional  plot 
that  involves  the  four  variables.  The  most  straight  forward 
approach  is  to  define  the  50*  probability  for  the  mean  of 
the  means  (Figs.  21  and  22)  in  terms  of  a  surface.  A  plane 
would  be  the  simplest  surface  choice.  The  equation  of  a 
line  perpendicular  to  the  50%  probability  plane  intersecting 
its  midpoint  in  the  distribution  planes  shown  would  then  be 
calculated.  All  of  the  precipitation  and  no-precipitation 
points  would  be  projected  onto  the  line  and  the  line  would 
become  the  x-axis  in  a  new  plot.  Thus  this  x-axis  reflects 
the  visual  and  infrared  mean  of  the  means.  The  visual  and 
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infrared  mean  of  the  standard  deviations  would  define  the 
y-axis  and  z-axis,  respectively. 

The  defined  plot  for  the  precipitation  and  no¬ 
precipitation  data  would  show  the  relative  dependence  of  the 
data  on  the  four  variables.  if  the  plots  produced  distinct 
groups  for  the  precipitation  and  no-precipitation  data,  new 
precipitation  probabilities  would  be  calculated.  The  varia¬ 
tion  with  satellite  resolution  size  of  the  precipitation 
probabilities  would  then  have  to  be  reconsidered  with  the 
new  data.  The  data  processing  described  is  beyond  the  scope 
of  this  particular  research  effort  and  is  reco mm mended  for 
future  investigation. 

Precipitation  specification  for  overcast  ceil¬ 
ings  can  be  delineated  by  the  values  for  the  naan  of  the 
means  and  the  mean  of  the  standard  deviations  for  any  array 
size  discussed  in  this  study.  The  mean  of  the  mean  values 
when  used  with  Figs.  21  and  22  will  indicate  the  probability 
of  precipitation,  given  a  similiar  time  of  year  and  climato¬ 
logical  area.  A  simple  linear  bi-spectral  threshold,  based 
on  a  50%  probability  of  precipitation,  is  defined  approxi¬ 
mately  as  extending  from  a  1.00  albedo  and  a  290K  cloud  top 
temperature  to  a  0.60  albedo  and  a  210K  cloud  top 
temperature. 
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2.  Overcast  and  Broken  Ceil incs 

Once  again  the  question  arises,  are  the  precipita¬ 
tion  overcast  and  broken  (category  IB)  and  no-precipitation 
overcast  and  broken  (category  2B)  data  sets  sufficiently 
separated  to  allow  differentiation  of  the  two  populations? 
If  so,  how  much  overlap  is  there  between  the  two  data  sets? 
a.  bean  of  the  Means 

The  mein  of  the  means  statistics  for  the  precip¬ 
itation  overcast  and  broken  (category  IB)  versus  no-precipi¬ 
tation  overcast  and  broken  (category  2B)  data.  Table  XIII, 
show  there  is  a  .309  and  27. OK  difference  between  the  two 
populations  for  the  10  x  10  pixel  size  and  a  .339  and  29. IK 
difference  for  the  4x4  pixel  size.  The  respective  differ¬ 
ences  are  approximately  equal  to  one  and  one-half  standard 
deviations  of  the  means  of  either  of  the  two  populations. 

For  these  data,  the  differences  in  the  trends  cf 
the  mean  of  the  means  between  the  10  x  10  and  4  x  4  array 
sizes  for  the  precipitation  overcast  and  broken  reports  are 
more  dramatic  than  for  the  precipitation  overcast  reports. 
There  is  a  0.036  visual  and  a  2.2K  infrared  difference  in 
the  precipitation  overcast  and  broken  data  and  only  a  0.006 
visual  and  a  0.1K  in  the  no- precipxtat ion  overcast  and 


TABLE  XIII 

Precipitation  Specification  Overcast  and  Broken  Ceilings 


Precipi¬ 

tation 

No -Preci¬ 
pitation 

Precipi¬ 

tation 

No-Preci¬ 

pitation 

10  x  10 

1.0  x  ig 

4x4 

4X4 

Mean  of 
Means 

(VIS) 

.  580 

.  271 

.616 

.  277 

(IB) 

252.  8K 
(-20°C) 

27  9.8K 
(♦  7°C) 

250. 6K 
( -23°C) 

279. 7K 
(+7°C) 

St  andard 

De  viations 
of  Means 

(VIS) 

(IB) 

.  210 

20.5K 

.  189 

1  6.6K 

.225 

21. 5K 

.207 

17. 8K 

Mean  of  (VIS) 

.  173 

.  128 

.124 

.  094 

Standard 

Deviations  (IR) 

10.  4  K 

5.9K 

5.6K 

3. 2K 

Standard 

Deviation  (VIS) 

.  078 

.  064 

.072 

.  057 

of  the 

St  andard 

Deviations  (IR) 

7.  OK 

5.7K 

4 . 9K 

3.7K 

broken  data  between  the 

10 

x  10  and  4 

x  4  array 

sizes. 

Therefore,  variation  of 

the 

re  solution 

size  and  the  trend 

for  the  mean  of  the 

means 

can 

be  used. 

as  well  as 

the  mean 

of  the  means  value 

itself 

»  in 

de lineating 

precipitation  from 

no-precipitation  for  these 

overcast  and  broken 

ceiling 

reports. 

b.  Mean  of  the  Standard  Deviations 


The  values  of  the  mean  of  the  standard  devia¬ 
tions  are  important.  The  0.045  albedo  and  4.5K  cloud  top 


temperature  difference  between  the  two  data  sets  for  the  10 
x  10  array  size  are  significant  as  well  as  the  0.030  and 
2.4K  difference  for  the  4x4  array  size. 

c.  Standard  Deviations  of  the  Means 

The  differences  in  the  standard  deviations  of 
the  means  are  nearly  equal  when  comparing  the  two  10  x  10 
array  sizes  and  the  two  4x4  array  sizes.  The  visual  dif¬ 
ferences  are  0.021  for  the  10  x  10  ana  0.018  for  the  4x4 
array  size.  Similiarly,  the  infrared  differences  are  3.9K 
for  the  10  x  10  and  3.7K  for  the  4x4  array  size.  The  vis¬ 
ual  differences  are  40%  less  than  the  statistically  signifi¬ 
cant  visual  mean  of  the  standard  deviation  differences  and 
the  infrared  differences  are  comparable  to  the  infrared  mean 
of  the  standard  deviation  differences. 

d.  Standard  Deviation  of  the  Standard  Deviations 


The  differences  in  the  standard  deviations  of 
the  standard  deviations  are  approximately  equal.  The  visual 
differences  are  0.0.14  for  the  10  x  10  and  0.015  for  the  4  x 


4  array  size.  The  infrared  differences  are  1.3K  for  the  10 
x  10  and  1.2K  for  the  4x4  array  size.  Both  rhe  visual  and 
infrared  differences  are  50%  less  than  the  respective  sig¬ 
nificant  differences  in  the  mean  of  the  standard  deviations. 
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Figure  26.  No-Precipitat icn  Overcast  and  Broken  for  4x4 
Array  (The  mean#  .277  and  279. 7K,  interval  is 
boxea.  The  2%.  3%.  5%,  7%,  and  \2%  frequencies 
are  for  147,  221,  j68,  515,  and  883  occurrences, 
respectively. ) 

f.  Precipitation  Probabilities 

Precipitation  probabilities  (Figs.  27  and  28) 
were  computed  from  the  precipitation  overcast  and  broken 
(category  IB)  and  the  no-precipitation  overcast  and  broken 
(category  2B)  cata  for  the  10  x  10  (Figs.  23  and  24)  and  the 
4x4  (Figs.  25  and  26)  array  sizes.  As  in  the  overcast 
ceiling  reports,  there  is  a  greater  than  50%  probability  of 
precipitation  at  the  low  albedo  values,  0.00-0.  20.  In  Fig. 
27,  the  75%  probability  between  230K-250K  interval  results 
from  three  precipitation  reports  of  four  total  reports.  The 


overcast  ceiling  reports  account  for  two  of  the 
precipitation  reports.  Therefore,  one  precipitation  and  one 
no-precipitaticn  broken  report  have  been  added  to  the  inter¬ 
val.  Similiarly  in  Pig.  28,  the  67%  probability  in  the 
210K-230K  interval  results  from  two  precipitation  reports  of 
three  total  reports.  The  overcast  ceiling  reports  account 
for  two  of  the  precipitation  reports.  Thus,  one  no-precipi¬ 
tation  report  has  been  added  to  the  interval.  A  few  reports 
in  these  low  albedo  and  cold  cloud  top  temperature  intervals 
are  producing  misleadingly  high  precipitation  probabilities. 
The  50%  probability  line  is  nearly  constant,  at  a  0.80 
albedo,  at  all  cloud  top  temperatures  for  the  overcast  and 
broken  ceiling  data. 

The  same  treatment  for  displaying  the  four  vari¬ 
ables  in  a  three  dimensional  diagram  suggested  in  the  previ¬ 
ous  section  (IV. C.  1.),  is  recommended  for  these  data  as 
well.  The  new  displays  then  could  be  used  to  calculate  pre¬ 
cipitation  probabilities  that  might  allow  a  more  accurate 
bi-spectral  threshold  specification  of  precipitation, 
g.  Summmary 

As  for  the  overcast  ceilings,  precipitation 
specification  of  the  overcast  and  broken  ceilings  can  be 
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in  this  study.  The  mean  of  the  means  values  when  used  with 
Pigs.  27  and  28  will  indicate  the  probability  of  precipita¬ 
tion,  given  a  similiar  time  of  year  and  climatological  area. 
&  very  simple  threshold  fcr  these  data  is  actually  dependent 
upon  the  visual  data  value  being  greater  than  0.80. 

D.  CONVECTIVE  VERSUS  CONTINUOUS  PRECIPITATION 

The  physical  processes  involved  in  convective  precipita¬ 
tion,  or  shcwers,  are  different  from  those  processes  usually 
involved  in  continuous  precipitation.  The  inherent  differ¬ 
ences  in  the  processes  might  result  in  a  statistical  separa¬ 
tion  in  the  thresholds  between  the  two  types  of 
precipitation.  For  this  data  set,  426  of  the  538  cases  are 
classified  as  convective  by  the  surface  weather  report. 
With  the  data  set  consisting  of  southeastern  United  States 
stations  in  August,  the  dominance  of  the  precipitation  cases 
by  convective  reports  is  not  surprising. 

1  •  SfiSfi  Statistics 

The  statistics  for  the  convective  and  continuous 
precipitation  (Table  XIV)  indicate  differences  in  the  means 
of  the  cell  means  of  .013  and  5.7K  for  the  10  x  10  array 
size  and  differences  of  .015  and  6.8K  for  the  4x4  array 
size.  As  mentioned  in  the  life  history  methods,  convective 
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precipitation  is  expected  to  be  associated  with  relatively 
higher  albedos  and  colder  cloud  top  temperatures.  For  these 
data,  the  convective  precipitation  are  associated  with 
colder  cloud  top  temperatures  resulting  from  the  cumulonim¬ 
bus  clouds. 

TABLE  XIV 

Continuous  versus  Convective  Precipitation  Specification 


Contin¬ 

Convec¬ 

Contin¬ 

Convec¬ 

uous 

tive 

uous 

tive 

12.  £  12 

10  x  10 

4x4 

4x4 

—  — 

Mean  of  (VIS) 

Means 

.  58  9 

.  576 

.602 

.617 

(IH) 

257.5K 

25  1.8K 

256. 2K 

249. 4K 

(-  16°C) 

(-2  1°C) 

(-  17°C) 

(-24  °C) 

Standard  (VIS) 
De  viations 

.  201 

.  214 

.221 

.228 

of  Means  (IR) 

1  9.  4  K 

20. 8K 

20.  IK 

21 .9K 

Mean  of  (VIS) 

Standard 

.  147 

.  180 

.099 

.  131 

Deviations  (IR) 

Standard 

7.6K 

1  1.  IK 

3.8K 

6.  OK 

Deviation  (VIS) 
of  the 

Standard 
Deviations  (IR) 

.075 

.  078 

.058 

.075 

5.3K 

7.2K 

3.  2K 

5.  2K 

The  sta 

tistics  for 

the  convective  and 

continuous 

precipitation  (Table  XIV) 

suggest  that 

for  visual 

satellite 

values,  the  mean  of  the  standard  deviations  is  the  best  sta¬ 


tistic  for  differentiation  of  these  two  data  types,  while 
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for  infrared  values,  the  Dean  of  the  means  is  the  best 
statistic.  The  use  of  two  different  statistical  types  to 
qualitatively  specify  continuous  precipitation  from  convec¬ 
tive  precipitation  is  unique  to  this  pair  of  convective  ver¬ 
sus  continuous  precipitation  data  comparison. 

The  infrared  mean  of  the  cell  means  show  that  the 
convective  precipitation  has  colder  cloud  top  temperatures 
due  to  the  greater  vertical  development  of  the  clouds.  The 
greater  vertical  variation  of  the  convective  precipitation 
clouds  is  seen  in  the  visual  '•texture'*  or  visual  mean  of  the 
standard  deviation  statistics.  This  "texture"  difference 
also  appears  in  the  infrared  values,  but  it  is  not  as  sig¬ 
nificant  as  the  differences  in  the  infrared  mean  of  the 
means.  The  lack  of  a  large  difference  in  the  visual  mean  of 
the  means  between  the  convective  and  continuous  precipita¬ 
tion  reflects  the  averaging  of  open  areas  and  cumulus  clouds 
associated  with  convective  precipitation.  Note  that  for  the 
fine  resolution  4x4  array  size,  the  visual  mean  of  the 
means  is  larger  for  the  convective  than  the  continuous 
precipitation. 

The  visual  means  of  the  standard  deviations  have  a 
difference  of  0.033  for  the  10  x  10  and  0.032  for  the  4x4 


array  size.  Contrast  these  difference  to  the  visual  mean  of 
the  mean  differences  of  0.013  and  0.015  for  the  10  x  10  and 


4x4  array  size,  respectively.  The  infrared  mean  of  the 
■eans  have  a  difference  of  5.7K  for  the  10  x  10  and  6.8K  for 
the  4x4  array  size.  These  differences  are  larger  than  the 
infrared  mean  of  the  standard  deviations  differences  of  3.5K 
for  the  10  x  10  and  2.2K  for  the  4x4  array  size. 

2.  Standard  Deviation  Stati  sties 

The  standard  deviation  of  the  means  and  the  standard 
deviation  of  the  standard  deviations  display  quite  similiar 
and  relatively  insignificant  differences.  The  visual  stan¬ 
dard  deviation  of  the  means  vary  by  0.013  in  the  10  x  10  and 
0.007  in  the  4x4  with  the  infrared  values  varying  by  1.4K 
for  the  10  x  10  and  1.8K  for  the  4x4  array  sizes.  The 
visual  and  infrared  10  x  10  and  4x4  array  size  differences 
in  the  standard  deviations  of  the  standard  deviations  are 
0.  003,  0.  017,  1.9K,  and  2. OK,  respectively. 

3 •  Distribution  Discussion 

The  distributions  of  the  convective  and  continuous 
precipitation  for  the  10  x  10  array  size  (Pigs.  29  and  30) 
and  for  the  4x4  array  size  (Figs.  31  and  32)  demonstrate 
the  similiarity  of  the  distributions.  This  similiarity 


disallows  quantitative  separation  of  the  convective  and  con¬ 
tinuous  precipitation. 


Figure  29.  Convective  Precipitation  Data  for  10  x  10  Array 
Size  (The  mean,  .  576  and  251.8K,  interval  is 
boxed.  The  255,  3%,  and  5*  frequencies  are  for 
9,  13,  and  21  occurrences,  respectively.) 


4 .  Summary 


Defining  the  qualitative  delineation  of  convective 


versus  continuous  precipitation  in  terms  of  the  shift  toward 


colder  cloud  top  temperatures  for  the  convective  precipita¬ 


tion  cases  is  one  available  delimiter.  The  second  qualita 


tive  separation  is  in  terms  of  the  shift  toward  higher 


albedo  values  of  the  mean  of  the  standard  deviations  for  the 
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CLOUD  TOP  TEMPEBATUUE  (K) 

Figure  32.  Continuous  Precipitation  Data  for  4x4  Array 
Size  (The  mean,  .602  and  256. 2K,  interval  is 
boxed.  The  2%,  3%,  and  53S  frequencies  are  for 
2,  3,  and  6  occurrences,  respectively.) 

convective  precipitation  cases  as  compared  to  the  continous 
precipitation  cases. 

The  convective  precipitation  cloud  top  temperature 
mean  of  251, 8K  (-21°C)  compares  quite  well  with  those  asso¬ 
ciated  with  life  history  methods  (Table  VII)  .  The  continu¬ 
ous  precipitation  cloud  top  temperature  of  257.  5K  (-16°c) 
compares  well  with  the  -16°C  threshold  of  Wylie  (1979)  (a 
life  history  threshold  for  Montreal,  Canada,  in  summertime) 
and  with  the  -12°C  bi-spectral  threshold  of  Muench  and  Kee¬ 
gan  (1979).  These  research  efforts  focused  on  mid-latitude 
precipitation  associated  with  fronts. 


E.  INTENSITY 

As  there  is  little  statistical  difference  between  the 
convective  and  continuous  precipitation  data,  the  gualita- 

V 

tive  analysis  of  precipitation  can  be  based  strictly  on  all 
light  or  all  moderate/heavy  precipitation.  The  light  pre¬ 
cipitation  cases  are  expected  to  have  a  lower  mean  albedo 
and  warmer  cloud  top  temperature  when  compared  to  the  moder¬ 
ate/heavy  precipitation  cases. 

1  •  Statistical  Discussion 

The  statistics  fcr  the  two  precipitation  intensi¬ 
ties,  Table  XV,  show  0.100  and  8.6K  differences  for  the  10  x 
10  array  size  and  .117  and  8. IK  differences  for  the  4  x  4 
array  size  in  the  means  of  the  means.  The  differences  are 
approximately  equal  to  one-half  a  standard  deviation,  an 
insufficient  separation  for  specification  of  these  two  pre¬ 
cipitation  classifications.  The  remaining  three  sets  of 
statistics  (Table  XV)  show  nearly  identical  values  in  the 
light  and  moderate/heavy  precipitation  for  the  two  array 
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TABLE  XV 

Light  versus  Hoderate/Reavy  Precipitation  Specification 


Light 

Bod era te/ 

Light 

floderate/ 

Heavy 

Heavy 

11  £  11 

11  i  11 

4  J  4 

tH 

Bean  of 
Beans 

(VIS) 

.565 

.665 

.598 

.715 

(IB) 

(- I^C) 

24  5.6K 
(-2  8*0 

251. 9K 
(-21*0 

243. 8K 
(-29*C) 

Standard 

Deviations 

(VIS) 

.  210 

.  201 

.226 

.204 

21. 3K 

of  Beans 

(IB) 

20.5K 

20. 4K 

21. 6K 

Bean  of . 

(VIS) 

.  172 

.  177 

.124 

.129 

Standard 

De  viations  (IS) 

10.  2K 

1  1.9K 

5.5K 

6. IK 

St  andard 

De  viation 
of  the 
Standard 

(VIS) 

.078 

.  081 

.072 

.079 

Deviations  (IR) 

7.  OK 

7. OK 

4.9K 

5.2K 

2.  Distribution  Discussion 

The  distributions  of  the  light  and  soderate/heavy 
precipitation  for  the  10  x  10  (Figs.  33  and  34)  and  the  4  x 
4  array  size  (Figs.  35  and  36)  allow  visual  confirmation  of 
the  overlap  between  the  two  intensities.  Basically  with 
only  74  cases  cf  eo derate /heavy  precipitation,  the  frequency 
contours  (Figs.  34  and  36)  reflect  a  noisy  distribution 
created  by  a  sorting  of  just  two  or  three  reports  into  a 
given  interval.  Additional  intensity  data,  specifically 
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aoderate/heavy  precipitation  cases. 


are  needed  before  a 


threshold  delineating  qualitative  intensities  can  be 
proposed. 


Figure  33. 


Light  Precipitation  Data^for^lO  x  10  Array  Size 
[T  be  sec* 


(The  Man,  .565  and  254. 2K,  interval  is  boxed. 
The  2%,  31,  and  5%  frequencies  are  for  9,  14, 
and  23  occurrences,  respectively.) 


For  the  10  x  10  and  4  x  4  array  sizes  aoderate/heavy 
precipitation  (Figs.  34  and  36),  the  area  encompassed  by  the 
7%  frequency  isopleth  at  the  higher  albedos  and  colder  cloud 
top  teaperatures  gives  an  indication  of  a  possible 
separation  between  the  two  intensities.  However,  the  biso- 
dal  nature  of  the  10  x  10  distribution  (Fig.  34)  and  the 
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three  relative  aarLaa  of  the  4  a  4  distribution  (Fig.  36) 
suggest  that  there  will  ba  a  significant  aaount  of  overlap 
between  these  two  qualitative  intensities.  Because  the  sta~ 
tistics  for  the  two  rainfall  rates  are  nearly  identical 
(Table  XT)  t  with  the  exception  of  the  aean  of  the  aeans,  the 
aethod  discussed  in  the  precipitation  specification  section 
(IT.C.)  which  involved  the  aean  of  the  standard  deviations 
to  aid  in  further  differentiation  of  two  data  types  cannot 
be  applied  to  these  two  intensities.  The  delineation  will 
have  to  be  derived  froa  the  aean  of  the  aean  statistics  and 
their  resultant  distributions. 
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3.  SJllAIi 

There  are  slight  Manifestations  of  the  expected 
downward  shift  in  the  albedos  and  the  earner  cloud  top  temp¬ 
eratures  for  the  light  precipitation  relative  to  the  moder- 
at  e/heavy  precipitation.  However,  additional  data  are 
reguired  before  a  threshold  can  be  proposed  to  delineate 
light  versus  no  derate/heavy  precipitation.  The 

investigation  of  the  additional  data  should  begin  with 
consideration  of  the  aean  of  the  Beans  and  their 
distributions. 

During  the  cell  resolution  discussion  of  section 
XV.B  •  1. ,  a  significant  nuaber  of  low  albedo  values  (less 
than  0.40  albedo)  were  noted.  k  comparison  of  the  light 
precipitation  data  (Fig.  33)  and  precipitation  data  (Fig. 
10)  for  the  10  x  10  array  size  reveals  that  96%  (105  of  109 
reports)  are  classified  as  light  precipitation.  Further 
analysis  shows  that  79%  (86  of  109  cases)  are  classified  as 
light  convective  precipitation.  (The  decomposition  of  con¬ 
vective  and  continuous  precipitation  data  into  light  and 
moderate/heavy  are  not  shown.)  These  low  albedo  values  are 
not  unreasonable  for  light  convective  precipitation  which 
can  be  produced  by  small,  isolated  cumulus  congestus  clouds. 


wt-w  a  m  wwy  >ti  71 


V.  SOHHABT  AMD  COHCLOSIOMS 

A.  DATA  PB0C2S5IMS  SOMHABT 

The  data  were  first  filtered  for  the  correct  time, 
1200-2000  GMT,  to  ensure  reliable  visual  satellite  values. 
The  visual  counts  were  converted  to  albedos  according  to  the 
duench  and  Keegan  (1979)  noraalization  scheae  (see  Appendix 
A)  .  This  scheae  corrects  for  Laabertian  scattering  as  well 
as  anisotropic  cloud  radiation.  Sorting  the  data  according 
to  the  Service-A  reports  (current  weather  and  cloud  group) 
provided  seven  precipitation  groups:  convective,  continuous, 
light,  aoderate/heavy,  general,  overcast  ceiling,  and  over¬ 
cast  and  broken  ceiling  precipitation  and  two  no-precipita- 
tion  groups:  no-precipitation  overcast  and  no-precipitation 
overcast  and  broken.  In  order  to  investigate  the  iapact  of 
satellite  resolution  on  precipitation  specification,  the 
stated  groups  were  subdivided  into  four  sizes:  10  x  10  (45  x 
45  nai)  ,  8  x  8  (36  x  36  nai),  6  x  6  (27  x  27  nai) ,  and  4x4 
(22  x  22  nai)  pixel  array  sizes.  The  aeans  of  the  albedos 
and  cloud  top  teaperatur es  for  each  array  size  for  each 
weather  or  cloud  group  were  calculated.  The  Bean  values 
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war®  sorted  into  10K  cloud  top  temperature  intervals  and 
0.10  estisated  albedo  intervals  to  provide  a  visual  descrip¬ 
tion  of  the  distributions.  The  standard  deviation  of  the 
■eans  and  the  aean  and  standard  deviation  of  the  standard 
deviations  for  both  albedo  and  cloud  top  temperature  were 
also  coaputed  for  each  distribution. 

B.  STATISTICS  SUHMARY 

The  coeputed  aean  and  standard  deviation  of  the  aeans 
and  the  aean  and  standard  deviation  of  the  standard  devia¬ 
tions  for  both  albedo  and  cloud  top  teaperature  were  evalu¬ 
ated  for  specification  of  precipitation,  specification  of 
continuous  versus  convective  precipitation,  and  specifica¬ 
tion  of  light  versus  so  derate  /heavy  precipitation. 

For  precipitation  specification  (Tables  III  and  XIII), 
the  visual  and  infrared  aean  of  the  aeans  produced  the  larg¬ 
est  statistical  differences  between  the  precipitation  and 
no-precipition  data  sets  for  both  the  overcast  and  the  over¬ 
cast  and  broken  ceiling  reports.  For  the  overcast  ceiling 
reports  (Table  XII),  the  separation  in  the  visual  and  infra¬ 
red  aeans  of  the  aeans  were  one  standard  deviation.  For  the 
overcast  and  broken  ceiling  reports  (Table  XIII) ,  the  sepa¬ 
ration  in  the  visual  and  infrared  aeans  of  the  aeans  were 


on*  and  one-half  standard  deviations.  The  aean  of  the  stan¬ 
dard  deviations  and  the  trend  of  the  aean  of  the  Beans  also 
yielded  a  large  difference  between  the  precipitation  and 
no-precipitation  data  sets. 

For  specification  of  continuous  versus  convective  pre¬ 
cipitation  (Table  XIV),  the  infrared  aean  of  the  aeans  and 
the  visual  aean  of  the  standard  deviations  yielded  the  larg¬ 
est  statistical  differences.  The  convective  precipitation 
infrared  seen  of  the  aeans  was  6K-7K  colder  than  the  contin¬ 
uous  precipitation  value  and  the  convective  visual  aean  of 
the  standard  deviations  was  20%  greater  than  the  continuous 
precipitation  value.  These  two  statistics,  however,  were 
only  indications  of  possible  separation  as  the  two  data  dis¬ 
tributions  displayed  significant  overlap.. 

For  light  versus  aoderate/heavy  precipitation  specifica¬ 
tion  (Table  XV),  the  aean  of  the  aeans  yielded  the  largest 
statistical  differences.  However,  the  relatively  saall 
quantitative  differences  (one-half  standard  deviation)  in 
the  visual  and  infrared  values  only  provided  an  indication 
of  the  upward  shift  in  albedo  and  the  downward  shift  in 
cloud  top  teaperature  for  the  aoderate/heavy  precipitation 
relative  to  the  light  precipitation. 
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C.  PRECIPITATION  PROBABILITY  StJ HHARY 


For  the  overcast  calling  reports,  the  10  x  10  array  size 
(Fig.  21)  50%  probability  line  decreased  linearly  to  lower 

albedos  with  colder  cloud  top  tesperatures ,  while  the  4x4 
array  size  (Fig.  22  )  50%  probability  line  decreased  linearly 
to  approximately  240K  and  then  increased  linearly  at  colder 
cloud  top  teeperatures.  For  the  overcast  and  broken  ceiling 
reports,  the  10  x  10  array  size  (Fig.  27)  50%  probability 

line  was  constant  at  0.80  albedo  while  the  4x4  array  size 
(Fig.  28)  50%  probability  line  gradually  decreased  linearly 
to  approximately  240K  and  then  gradually  increased  linearly 
at  cclder  cloud  top  temperatures. 

0.  DATA  DISTBIB0TI0H  SO HHARY 

The  distributions  of  all  seven  precipitation  classes  and 
both  no- precipitation  classes  were  non-Gaussian.  This 
result  disagrees  with  the  assertion  of  Lovejoy  and  Austin 
(1979)  that  their  cumulus  rain,  non-cumulus  rain,  and  non- 
cumulus  no-rain  data  distributions  were  Gaussian.  %  At  the 
relatively  fine  resolution  4x4  array  size,  all  seven  pre¬ 
cipitation  classes  and  the  no-precipitation  overcast  data 
displayed  bimodal  distributions  while  at  the  relatively 
coarse  10  x  10  array  size,  all  displayed  an  uninodal 
distribution. 
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E.  CONCLUSIONS 


The  conclusions  f roe  this  study  of  precipitation  speci¬ 
fication  are: 


Varying  the  satellite  data  resolution  fros  484  nmi*  (4 
z  4  array  size)  to  2025  nmi*  (10  x  10  array  size) 
results  in  a  statistically  significant  difference  in 
the  representation  of  precipitation  or  no-precipitation 
data.  Variation  in  the  distribution  functions  and  the 
characteristic  aeans  and  standard  deviations  with 
increasing  cell  size  deaonstrates  systematic  trends 
which  may.  with  further  study,  provide  improved  basis 
for  rain  detection  and/or  quantification. 


tion  is  defined  as  extending  from  an  albedo  of  1.00  and 
a  cloud  top  teaperature  of  290K  to  an  albedo  of  0.60 
and  a  cloud  top  teaperature  of  210K  for  overcast  ceil¬ 
ing  reports.  For  overcast  and  broken  ceiling  reports, 
a  threshold  based  on  the  albedo  being  greater  than  0.80 
specifies  a  50%  probability  of  precipitation.  However, 


specifies  a  50%  proba 
the  precipitation  pro 
lite  resolution  size. 


a  50%  probability  of  precipitation.  However, 
itation  probabilities  varied  with  the  satel- 


•  The  precipitation  and  no-precipitation  data  sets  can  be 
differentiated  by  the  mean  of  the  aeans,  the  mean  of 
the  standard  deviations,  and  the  trend  of  the  mean  of 
the  aeans. 


•  Differentiation  of  convective  precipitation  from  con¬ 
tinuous  precipitation  shows  proaise  through  the  differ¬ 
ences  in  the  infrared  aean  of  the  aeans  and  through  the 
visual  aean  of  the  standard  deviations. 


•  Qualitative  specification  of  light  versus  moderate/ 
heavy  precipitation  shows  soae  promise  through  a 
threshold  based  on  the  relative  shift  toward  colder 
cloud  top  temperatures  and  higher  albedos  of  the  noder¬ 
ate/heavy  precipitation  relative  to  the  light  precipi¬ 
tation. 


F.  SUGGESTED  FURTHER  STUDI 


The  recommendations  for  further  study  are: 


•  The  proposed  data  processing  that  would  display  the 
visual  and  infrared  aean  and  standard  deviation  of  the 
means  (ZV.C.1.)  should  be  investigated  for  both  the 
overcast  and  the  overcast  and  broken  data  to  determine 
if  the  precipitation  and  no-precipitation  data  can  be 
more  distinctly  separated. 


•  The  appearance  of  biaodal  distributions  at  the  fine 
resolutions  (6x6  and  4x4  array  sizes)  in  the  seven 
precipitation  categories  and,  in  the  no-precipitation 
overcast  category  should  be  investigated  as  a  possible 
synoptic  signature. 

•  Poraal  discriminant  analysis  should  be  applied  to  the 
relevant  statistics  to  yield  confidence  levels  of  indi¬ 
vidual  visual  and  infrared  satellite  data  pairs  for 
precipitation  specification. 


•  The  data  should  be  further  investigated  at  the  2  x  2 
array  size  to  analyze  this  resolution  m  representing 
the  data. 
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APPENDIX  k 


The  Huench  and  Keegan  (1979)  normalization  relates  the  nor¬ 
malized  reflectivity,  r”n  ,  to  the  varying  solar  angle  and 
maximum  digital  counts  through  the  reflectance  term,  'r',  and 
the  anisotropic  scattering  through  the  X  term.  Table  XVI 
defines  the  symbols.  Table  XVII  lists  the  geometric  identity 
equations  and  Table  XVIII  list  the  normalization  equations. 
Fig.  37  gives  an  example  of  the  normalization  applied  to  the 


stated  location 


TABLE  XVI 
List  of  Syabols 


Description 


Units 


GOES  video  count  number  (0-63) 

dimensionless 

GOES  video  count  number  for  perfect 
diffuse  reflector  and  overhead  sun 

dimensionless 

Greenwich  meridian  time 

hours-minutes 

seconds 

Distance  of  earth  to  sun 

km 

Mean  distance  of  earth  to  sun 

km 

Julian  date 

dimensionless 

Hour  angle 

radians 

Cloud  reflectivity 

dimensionless 

Arc-length  observer  to  subsatellite 
point 

radians 

Declination  of  the  sun 

radians 

Zenith  angle  of  the  sun 

radians 

Longitude 

radians 

Longitude  of  subsatellite  point 

radians 

Anisotropic  scattering  coefficient 

dimensionless 

Latitude 

radians 

Azimuth  of  the  sun 

radians 

Azimuth  of  the  satellite 

radians 

if 


TABLE  XVII 

Basic  Geometric  Satellite- Earth  Relationships  * 


Declination: 

4  -  0.408  sin  [(0-81)  *  2*/ 365] 

Solar  distance  ratio: 

8/B0  »  1  -  0.167  cos  [ ( d-14 )  *  2*/ 365] 

Hour  angle: 

h  *  A  ♦  *  -  G( hours)  x  ir/12 
ire- length : 

cosy  ”  cos(Ab  -  A)  cos* 

Satellite  azimuth : 

sin  (*2  -  *)  -  sin  (Aa  -  A)/siny 
Solar  azimuth  angle: 

cose  "  sin*  sin*  +  cos*  cos$  cosh 
Solar  azimuth: 

sln*1  •  cosd  slnh/sinc 


ingles  In  radians 


TABLE  XVIII 

fluench  and  Keegan  normalization  Equations  * 


?  -  <£-)2  *  seec 

uo 


a  * 


-  1+2  -  +xl 


C.  -  com-  CCC  -  50)  *  1.8) 


C2  •  0.7  cos  CCC  -  22.5)  *  4)  *  (1  -  cose) 


C-  -  cos°  ((A+  -  70)  *  1.3) 


X  8  1.0  +  0.05  *  (1  +  cos(2*e>)  +  0.20  *  (Cj  ♦  C2)  *  Cg 

*n  -  tl.09  -  2*(1.09  -  F  *  X  *(R/H0)2)/<1  ♦  cos1/2  C)J 


Angles  in  degress 


& 


X 


'  X 
iV 


\  % 

\  *, 
\ 

\ 

V 

\ 

.  \ 
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Figure  37.  Serialized  Cloud  Reflectivity  as  a  Function  of 

rideg  Count  and  Tiae  of. Day  for  April  16  at  42°H 
and  74®W,  with  Calibration  Count  (cj  egual  to 
6|^  (Ezaaple  taken  froa  Huench  and°Keeganf 
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